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TIME DEPENDENT MONOENERGETIC NEUTRON TRANSPORT IN A 


FINITE SLAB WITH INFINITE REFLECTORS 

By 

Perry Allan Newman 
ABSTRACT 

The initial-value transport problem of monoenergetic neutrons 
migrating in a thin slab surrounded by infinitely thick reflectors 
is solved using the normal -mode expansion technique of Case. The 
results obtained indicate that the reflector may give rise to a 
branch-cut integral term typical of a semi -infinite medium while the 
central slab may contribute a summation over discrete residue terms. 
Exact expressions are obtained for these discrete time eigenvalues 
and numerical results are presented showing the behavior of real 
time eigenvalues as a function of the material properties of the 
slab and reflector. These eigenvalues are finite in number and all 
of them may disappear into the branch cut or continuum as the material 
properties are varied; such disappearing eigenvalues correspond to 
exponentially time -decaying modes. The two largest eigenvalues can be 
compared with critical dimensions of slabs and spheres and it is shown 
that the numerical values agree with criticality results of others. 

In the limit of purely absorbing reflectors or a bare slab, the present 
solution has the properties which have been previously reported by 
others who used the Lehner-Wing technique to solve corresponding 
problems . • 
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I. INTRODUCTION 


A few time -dependent, raonoenergetic neutron transport problems in 
plane geometry have been solved recently by applying Case ' s normal-mode 
expansion technique (refs. 7) to the equation which results when the 
time dependence has been removed by a Laplace transformation. In these 
time -dependent solutions, contributions due to various parts of the 
spectrum of the transport operator have been Indicated by suitably 
deforming the integration contour of the inverse Laplace transformation. 
This approach was used by Bowden (refs. 1, 4) for his analysis of time- 
dependent, one -speed neutron transport in a bare slab of finite thick- 
ness with isotropic scattering, a problem which held been treated 
extensively by Lehner and Wing (refs. l6, 17)- Another successful 
application of this technique was made by Kuscer and Zweifel (ref. l4) 
to the time-dependent, one-speed albedo problem for a semi-infinite 
medium which scatters neutrons isotropically. For these two applications 
the time -dependent solution involves only a single material medium. In 
each of these problems, construction of Case's normal-mode expansion in 
the transform plane depends upon the two material properties which 
characterize a single uniform medium with isotropic scattering: the 

total macroscopic cross section, or, and the average number of 
secondaries per collision, c. For time -dependent problems in which 
more than one medium is involved, the transform plane must be taken as 
the superposition of "single -medium" planes, one for each medium. The 
situation then for a problem in which material properties vary from 
point-to-point will be very complicated. 
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Mika (ref. 18) has studied one such problem: the initial value 

problem for raonoenergetic neutrons in a nonuniform slab surrounded by 
a vacuum. He used the same approach that Lehner and Wing (refs. 1 6 , 17) 
had used for a uniform bare slab problem and, as might be expected, his 
more general hypothesis results in fewer details. In particular, it 
appears that theorems concerning the reality and number of discrete 
time eigenvalues cannot be established. At the outset, Mika (ref. 18) 
indicates that such results would be used, in practice, for a system of 
uniform slabs. Even for these cases in which there are a limited number 
of different material media, he states that the most suitable means of 
calculating discrete time eigenvalues would seem to be the normal-mode 
expansion approach employed by Bowden (refs. 1, k ) . This is the 
approach used in this thesis to analyze a simple idealized two-media 
problem in which one would expect to have discrete time eigenvalues in 
Order to obtain some insight concerning their behavior as a function of 
material properties. Such an approach has been utilized for one two- 
media time -dependent problem by Erdmann (refs. 8, 9) who investigated 
the time decay of a plane isotropic burst of monoenergetic neutrons 
introduced at the interface of two dissimilar semi-infinite media which 
scatter isotropically. In his solution, contributions due to the 
continuous spectrum are different for the two media; apparently the 
continuous spectrum depends on x. There are no discrete eigenvalues 
in his problem. 

Lehner (ref. 15) has demonstrated that the continuous spectrum of 
the transport operator is very sensitive to the explicit formulation of 
a physical problem. He considered a slab of finite thickness 

1 / 
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surrounded by a pure absorber which had the sane total macroscopic cross 
section as the slab. He obtained the same point spectrum as that found 
for the bare slab (refs, 16, 17) but found the continuous spectrum to be 
only the imaginary axis Instead of the entire left-half plane. Very 
recently Hintz (ref. 10) has generalized Lehner’ s problem by allowing 
the pure absorber to have any cross section and found that when the two 
total macroscopic cross sections were different (o^ ^ Bee Fig, 1) 
the continuous spectrum is a strip parallel to the imaginary axis of 
width | < 7 ^ - a 2 1 sued that the point spectrum may be empty. He shows 
that his results reduce to those of Lehner (ref. 15) when * Op 
but does not indicate how the bare slab results of Lehner and Wing 
(refs. 16, 17) can be recovered. In the present problem, a finite slab 
is surrounded with a material which can scatter as well as absorb 
neutrons. Thus the bare slab and slab surrounded by pure absorbers are 
special cases and it is shown that the present solution has the proper 
behavior (refs. 10, 15, 1 6, 17) for these special cases. 

Consider a slab of material which scatters neutrons isotropically, 
extends from x = -a to x = a and is characterized by the nuclear 
properties erg 811,1 c 2* uniform slab is surrounded by uniform 

infinitely-thick reflectors of another material characterized by the 
properties and c^ (see Pig. 1). For a physically meaningful 
system, these reflectors should be nonmultiplying media since they 
extend to infinity. Therefore, we take e-^ < 1. For isotropic 
scattering of monoenergetic neutrons in a sourceless medium and plane 
geometry, the neutron angular flux, i|r(x, U; t), satisfies the 
equation (ref, 7) 








If 
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+ n + a(x)t(x,n,t) =. £&§« f\( x ,»',tw, 

( 1 . 1 ) 

where t is the neutron speed multiplied by the real time, x and n 
are shown on Figure 1 while <j(x) and c(x) are given by 


cr(x), c(x) 


\ a l> C 1 for 
|cr 2 , c 2 for 


| x | > a 
J x | < a. 


( 1 . 2 ) 


We seek the solution of this equation subject to the boundary conditions 


. i 1 ijr(±x,n,t) = 0 (1.5) 

I X I — » oo 

and the continuity conditions 

^(ta+,|i,t) * \|f(±a-,|i,t), (1.4) 

given the initial condition 

t(x,H,0) * f(x,n) (1.5) 

which we assume satisfies (1.5) and is extendable without poles or 
branch cuts in the finite |i-plane except perhaps for a discontinuity 
across the imaginary axis. When the material properties of the 
reflectors are taken to be those of a vacuum this problem reduces to 
that of Lehner and Wing (refs. l6, 17) while for a pure absorber it 
reduces to that considered by Lehner (ref. 15) or Hintz (ref, 10). 

The method of attack to be used in solving this problem is the 




following: 
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1* Remove the t -dependence with a Laplace transformation, 

2. Solve the transformed equation "by applying Case's technique. 

3* Determine the analytic properties of this transformed solution 
in some right-half s -plane. 

4-. Recover the t-dependence and simplify by suitably deforming the 
integration path of the inverse transformation. Previously cited re- 
sults (refs. 1, k, 8, 9; 1^) lead us to expect that the reflectors 
should contribute continuous-spectrum type terms typical of a semi- 
infinite medium while the central slab should give rise to some point - 
spectrum type terras and their corresponding discrete time eigenvalues. 

5* Calculate real discrete time eigenvalues as a function of 
material properties if and when they exist. 

This is the method which has been successfully employed by Bowden 
(refs. 1, k) f Kuscer and Zweifel (ref. l4), and Erdmann (refs. 8, 9 ) > 
we use many of their results in solving the present problem. In fact, 
our solution contains parts which resemble their solutions. Some 
preliminary results for the present problem were given in ref. 23. 


gwMraiiM Mw i 




II. TIME REMOVAL AND ELEMENTARY SOLUTIONS 


If we take the Laplace transformation of \Jr(x,p.,t) as 


t s 



\|r(x,n,t)dt. 


( 2 . 1 ) 


then the inverse transformation required to recover the t-dependence is 

>7+iw 
y-ico 


n py+iw 

t(x^,t) = x-r- / e st \|r s (x,n)ds, 

^nx j y. 


( 2 . 2 ) 


where y is to the right of all singularities and "branch cuts of 
ty s (x,|i) in the transform plane; that is, the s-plane. From previously 
cited work of others, it is expected that the path of integration in 
Eq. (2.2) can be deformed to indicate more precisely the character of 
vjr(x,}i,t). When we apply the transformation of Eq. (2.1) to Eq. (l.l), 
integrate by parts in the usual manner and make use of the initial 
condition (1.5), we obtain 


[x t s (x,n) + [a + cr(x)j i s (x,n) ■ J* \lr s (x,n')<4i’ + f(x,n). 


(2.5) 


Equations (l.J) and (1.4) become under this transformation 


and 


\|r s (±x,|i) = 0 


(2.4) 


t s (±a+,n) = i|r s (±a-,|i). 


(2‘5) 
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Before applying Case's technique to solve Eq. (2.5) subject to 
conditions (2.4) and (2-5), let us examine some properties of the 
transformed solution which follow directly from the governing equations. 
Bowden (refs. 1 , 4) introduced these ideas at a later step in his work, 
but here they aid in the construction of the solution. 

An arbitrary function of two variables, f(x,u), can be written as 
the sum of its even and odd parts, namely, f + (x,ji) and f_(x,|i). They 
are given, of course, by 

f+(x,n) = | [f(x,|i) ± f(-x,-|i)J (2.6) 

and have the property 

f ± (-x,-n) =±f+(x,n). (2.7) 

Since c(x) and cr(x) are even functions of x, we can easily show 
from Eq. (2.5) that the even and odd parts of i|r g (x,|jO obey the 
equation 

\i ^\|r s +(x,|a) + [s + cr(x)] ^ s± (x,n) = ^ s± (x,h' )<4*’ + f±(x,|i). 

( 2 . 8 ) 

The boundary conditions for \|r s± corresponding to Eqs. (2.4) and (2.5) 
are written as 


v , ^ te(x,|i) = 0 

X I — » CO 


(2.9) 




and 
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where the ± subscripts denote definite parity parts of a function 
(see Eqs. (2.6) and (2.7))* Equations (2.8) -(2.10) tell us the following: 

1. All solutions of the homogeneous equation associated with (2.8) 
can be made to have a definite parity. 

2. The boundary conditions preserve the parity. 

j, The definite parity parts of an initial distribution excite 
inhomogeneous solutions of corresponding definite parity. Therefore, we 
can separate this problem into two problems, one for \{r s+ , the other for 

and combine the results at any stage of the calculation. 

The functions f+(x,fi) and t s ±(x,|i) are broken up as 


jf-L+U,^), |x| > a 

f±(x,|i) = / 

|^P4- (XjN’) } | x | a , 


( 2 . 11 ) 


and 


t s ±(x,h) = 


Vl±(x,»i,s), |x| > a 

t2±( x ^h; s), | x | < a. 


( 2 . 12 ) 


so that Eqs. (2.8), (2.9), and (2.10) become 


8 

U ££ fj+(x,H,s) + (s + Oj)<|Tj ± (x,H,s) = J 


+ f J± (x,n), J = 1,2, (2.13) 


and 


|x| ^ ^± (x ^' s) ■ 0 


= t2±( a >H>s) 


( 2 . 11 ;) 


(2*15) 
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The notation gj+(a,n) means the limit of g£(x,n) as x -» a from 
r )diura j. Solutions of Eqs. (2.13) will he obtained by constructing 
even and odd particular solutions, tjp±(x,n), and adding to them 
solutions of the corresponding homogeneous equations, \Jrj c+ (x,|i), 
so that conditions (2.14) and (2.15) can be satisfied. These functions 
and will be constructed from Case's elementary solutions 

which we shall denote here as \|rj v (x,|i, s) . 

The elementary solutions, \|/\j v (x,|i,s), are solutions of the 
equation 

^ ■ jL P 

H ^t Jv (x,n,s) + (s + <yj)ij V (x,|i,s) * - CjOTj J tjv^fiSsJdii* 

“1 


in the form 


tj v (x,n,s) = 9j SV (^) e "^ s+a ^ x ^ v 


where v is a complex parameter introduced in this separation of 
variables and sv (m-) is normalized as 

f * s + <Jy 


Bowden (refs. 1, 4) and Erdmann (refs. 8, 9) have investigated these 
solutions; many of their results are given in Appendix A and will be 
used herein. They show that the solutions <Pj sv (u) are given by 

‘P.jsvM = I c J®j vP 7^7 + *j S ( V ) 8 ( V - l*)> V 6 (-1,+1), 


(2.19) 
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where P denotes the Cauchy principal value, 5(v - u) is the Dirac 


delta function and 


?^ s (v) * s + cfj - tanh” 1 v, 


and two discrete solutions. 


fevojM 


1 c j a ^) v 0j 

? v — rr? s € s ji' 

* v oj + * 


( 2 . 20 ) 


( 2 . 21 ) 


provided that the function flj s (z), 

nj s (z) = s + ffj - 


1 1 

Cj0jZ tanh — , 


(2.22) 


of two complex variables s and z vanishes at the two points ±v Q j . 
The condition for this to happen (refs. 1, 4) is that s lie inside 


the curve C* (s e S** , see Pig. 2) defined by 


Is + CT^ PS' 

C 1 = < v = a’ + ip’ a* = 

J \ c j ff j * 


tanh 


■W) • 


(2.23) 


We note that Vq^ is an analytic function of s for s e except 

for a branch cut on the real s-axis between -aj and -aj(l - cj). We 
have denoted by + v qj that zero of flj s (z) for which Re(vQj) >0, 
s 4 ^-aj,-(Tj(l - Cj)J . The important result is that the general solution 
of (2.16) can be expressed as the linear combination 

>j(x,|i,s) « |a^v 0 j(x,^s) + bjt-v 0 j( x ^s)^6j(s) 

+ f 1 Aj(v)\|rj v (x,|i,s)dv, (2.24) 



15 


where we define 6^(s) as 


Ms) 


- f 1 * 

X.0, 


s e S. 


s e Sj e 


and the s -dependence of the expansion coefficients has not been 
indicated. To avoid confusion with our notation ± for parity, the 
customary notation for the discrete modes has not been used. 

The solution (2.24) does not have definite parity. For a medium 
which is connected and symmetric about x m 0 (such as our slab), even 
and odd solutions can be written as 

tj+(x,d,s) = aj ± ^ Vq ^(x,h,s) ± ty-v 0 j s ) jS j ( s ) 


•f A < j ± (v)^ JV (x,n,s) ± tj(_ v )( x >^ s )J dv > 


where 


!± “ | [ a J * b j] 


j±C v ) * 


We have used the properties 


t±v 0 4 ( -*> -Vt s ) = %y 0 j ( x > s ) 


t J (+ v ) ( “ x > s) = (q:v) ( X >H> s ) 
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For a medium which extends to infinity in the x -direction (such as our 
reflector) the boundary conditions (2.4) require that 

bj m Aj(-v) S' 0, 0 < v < 1 if x -4*w 

or 

aj s Aj(v) sO, 0<V<1 if x-».« (2.29) 

for the expansion coefficients in Eq. (2.24) when Re(s) > -<jj. 

We use results (2.24), (2.26), and (2.29) now to construct the 
solutions tjc± 


III. CONSTRUCTION OF TRANSFORMED SOLUTION 


The even and odd homogeneous solutions in the slab, w can be 
written in the form of Eq. (2.26) with j » 2. On the other hand, the 
homogeneous solution in the reflectors, subject to the boundary 
condition (2.29)> can be witten from Eq. (2.24) as 

* f A 1 (v)t lv (x,n,s)dv, x < -a 


ti(x,H,s) * 


a-L’ty (x,(i,s)5 x (s ) +f A 1 ’(v)t lv (x,H,s)dv, x > a, 


(5-i) 

for Re(s) > -o^. The continuity conditions (2.5) and the parity of the 
solutions ^2c+ can te used to relate the coefficients in Eq. (3*1) • 

We find that an even solution inside the slab requires 

a l' = **1 and ( v ) ® A i(” v )> 0 < v < 1, (5* 2a] 


while an odd solution inside the slab requires 


a l' = -b^ and A^'(v) = -Ai(-v), 0 < V < 1. 


The explicit forms of ^2c+ an< ^ ^lc± are therefore 


to«+(x,^i,s) * a. 


2± [*v 02 <*^ 8 ) ± t.y 02 (x,(i,s)] 5g(s) 

f 0 k 2t^ ^2v^^ 8 ^ ± ^2(-v)( x '^ s j| 


(5-3) 


1 6 




i ' a a±^-v 01 ( x ^^ s ) 5 i( s ) 
±*Kf v IwW*) 


Hi 

+ A a +(-v)f 1 (. v )(x,n,s)dv, x < -£ 

± J A l+(-v)f lv (x,H,s)dv, x > a, 


for Re(s) > -o^. 

We turn now to construction of tjp+‘ Consider a function 
gj s (x,n;x 0 ) which satisfies the equation 


(3A) 


11 S’ g J s( 


*,m*o) + (s + 0 j)g Js (x,M.;x o ) 



gj S (x^V*o)<V 


+ 8(x - x 0 )fj(x 0 ,,u). (3-5) 


Upon integrating on x from Xq - € to Xq + € and taking the limit 
€ -» 0, we obtain the jump condition 


fj(x 0 ,fi) 

6j S (xo+,^;x 0 ) - gjstXQ-^xo) * — 


(3-6) 


The function ^^(x^ii^s) defined as 




(5- 7a) 


where 


tjp(x,li,s) » J gj S (x,ii)X 0 )dx 0j 


(3-7b) 


(medium) 

J 


IT 


is seen to be a solution of Eq,. (2.15)* It is shown in Appendix B that 
explicit forms of tyg>p4- and tip+ can te written as 

^2p±( x ^ ,s ) = rF S f(x,v 02 ,s)tv 02 (3c>^^ s ) 

± -x, v 02 , s)^_ VQ2 (x,|a, s)J 6g(s) 

pi 

+ J Q F 2f( x ^^s)f 2v (x,n,s) dv 

pi 

± J F 2± (-x,v, s)\|r 2 ( _ v ) (x,|i, s)dv, (5.8) 


tlp±( x >H> s ) = v oi> s )tv 01 ( x ^^ s ) 

+ ^ F i±( x ;-Voi^s) - F±(-a,v 01 ,s)J \|r_ V()1 (x^,s)| d^s) 
+ f Fi±(x, V, s)\Jr lv (x,|i, s)dv 

J 0 

+ / o [»*<*, -v ,s) - P ± (-a,v,s)J ti(„ v )(x,H,s)dv, 


x < -a, (3* 9a) 


for Re(s) > -cf^, and 


t 1]E ^(x,|x,s) 


= |-F±(-a,v 01 ,s) + F 1 ±(-x,-v 0 ps)j t V()1 (x,|a,s) 

± F 1 +(-x,v 01 ,s)t.v 0 i( x ?^ s )J 8 i( s ) 

± /. jjF+(-a,v,s) + P^C-x^-v^sij i|r lv (x,n,s)dv 


* J o Fq+( -x,v, s)ti( „ v ) (x,n, s)dv, x > a, (3-9b) 
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for Re(s) > where 


F ± (-a,ai,s) a P 1± (-a,-co,s) + F^-a^s), 

F2 ±(x, 0 ),s) a dx Q 

u -a 

F;i±(x,a>,s) a J C- L± (x 0; a))e^ s+<Jl ^ Xo ^ dx 0 . 


Here the C^ + are full-range expansion coefficients of the function 


fj+(x,n)Ai and are given "by 


: J ±(x 0’ v) ’ ■^TTTTTt/-, 

vflJ s (v)Oj S (v) J -1 


and if s e S^, 


C j±( x 0^ v O,j) “ 


c J a j v 0j 2 n , js (v 0 j) 


J_ x fj±(xo,t*)i> Voj (ii)ai* 


c i+(x0;“ v 0.l) = 


^r—Jl 


J M v OJ » js^” v Oj 


Throughout, we shall use + and - superscripts to denote the limiting 
values of a function on its "branch cut as the argument approaches the 
cut from the upper (+) and lower (-) half- planes. The function flj s (z) 
of Eq. (2.22) has a "branch cut along the real z-axis (-1,1) such that 


fl^g(v) = Aj S (v) ± irtCjffjv/2, -1 1 v < 1. 
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The functions n' j S (z) are defined by 

n 'js (z) 3 k a i^ z) - 


(3.13) 


Note that the parity of the coefficients Cj + is opposite that indicated 
by the ± subscript. Nevertheless, the solutions are easily seen 
to have the indicated parity. 

The solutions ijr^ (Eq. (2.12)) of our problem are written in 
terms of and as 


/ x| >a 

t s +(x^) « l 

Lt 2c +( x ^^ s ) + ^2p±( x >^ s )> | x I < a * 


(3.14) 


The solutions in medium 1 , Jx J > a, have been constructed so that the 
boundary condition (2«l4) is satisfied. Application of the continuity 
condition (2.13) allows us to determine the unknown expansion coefficients 
of which appear in Eq. (3-l4). That is to say, if we substitute 

x = a in Eq. (3*l4), apply the continuity condition (2.15) and use the 
explicit forms of given by Eqs. (3*3) and (3.4), we obtain a two- 

media full-range expansion involving the 9j SV which contains unknown 
coefficients aj+ and Aj + . The same expansion is, of course, obtained 
for x = -a. This type of expansion and its orthogonality relations are 
discussed in Appendix C and we show in Appendix D that such an expansion 
is obtained for the present problem. Erdmann (ref. 8) proved completeness 

vv 

theorems which apply in such time -dependent problems while Kuscer, 
McCormick and Summerfield (ref. 13) derived orthogonality relations 
which are applicable to two-media expansions which arise in time- 
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independent problems. In Appendix C, we extend their results to obtain 
orthogonality relations in a form which are valid for all regions of the 
transform plane. As usual in problems involving a slab, we cannot obtain 
closed form solutions for the expansion coefficients. However, we can 
use the orthogonality relations (App. C) to obtain expressions which 
give the expansion coefficients implicitly. That is, the continuum 
coefficients A2 +(v) are given as the solutions of Fredholm integral 
equations and all of the other coefficients are obtained from the 
AgfCv). More specifically, if we define 

E2±(v) ^ A 2t (v)fl£ s (v)S£ s (v)e( s+<T 2> a / v 

rnd 

E U (V) 8 A 1± (-v)^ 0 (v)£j£ s (v)e-( s+ ‘ , l) a / v , 


(5-15) 


the use of the orthogonality relations leads, after some algebra, to 
the following list of equations: 


Eg+Cv) = ^2±( v ) 


p+(v) ± f^22sC 
24 V 2 fi ls( 

fr 1 e-2( 

Vo -3 


k s Q 2s(°°) v / v s v\ 

2 % s (~) X ° ( , ] J 


-2(s+q 2 )a/p. Xo(-M,s)|i d^ 
n2sM&> s (v)(v + v) 


+ 82 (s)a^e-( s ^ a / V 0 2 Xo(- V0 2>s)— 

v 02 + V J 


(3.16) 


r- lICU 


J 


c 2 a 2 v 02^2s^ v 02^ a 2t^ S+a ^ a ^ 02 = J ^( v 02^ ± {if X 0^ V 02 ;S ^ 

f p 1 e“2(s-KT2)a/|i Xo(-|i^s)|i d|i 

X £0 E&t ^ «2s(» i ) fi 2s(^)(^ + v 02) 

+ |a^e-( s+ °2) B /v 02Xo( . Vo2)B) j, 


s e S, 


(5-17) 


E 1± (v) = I 3± (v) ± 2E fl | s(v)flis(v) E at (v)e- 2(s+ °2 )a/v 

^ fl| s (v)ni s (v) ** 


* (i ks 


f r 1 t \ e” 2 ( s+ff 2) a /M X 0 (-fi,s)cp lsv (n)2n dji 
x </_ — — — 


fi 2s^) J2 2s(^) c l ff l v 


+ 5g(s)ap4.e ( s+ar 2) a / v o2 Xq(-Vq2>s) 


V - V, 


0 < v £ 1 (5.18) 


and 
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1 

2 




x. - i *,<„ 4” aW “ '•<-»> * 



“ v 0l) 


+ 5 2 (s)a at e‘( s+a 2)a/vo2 X^-VQgjs) - 


v 02 


02 “ v 01l 


S € s^. 


( 3 . 19 ) 

The Ij + and J^ + terms contain only integrations over the initial 
distribution and are therefore known functions when f(x,n) is 
specified. They are given by 

W v > = ^Pa±(-a,v,s)e( s+<r l) a / v i£ s (v)Gg s (v) ± {| *s jjgpj- %( -v, s)J 

X < f f 1 P a (a,4,s)e-(«2>^ XqU.s) -BJt 

i y o u + v 

+ S 2 (s)T s± (e., V02 ,a)e- {s ' K2)&/v 02 Xc(-v 02 , B ) 


, «>u(-) f r 1 
+ JSPTL/o 


Fj+C-a^sJe 


(s+a 1 )a/n < l*2sv(l 1 ) 2 ^ ^4 
^©(^sjcgcrgv 


+ 6 1 (s)P it (.a,v 01 ,s)e( 6+ ®l) a / v 01 -—-i - 

AQV-V01,S) 


'01 


v - V 01 

0 ^ v ^ 1, (5.20) 
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J 2±^ v 02) 





%s(°o) 


F2fc(a,H,s)e-( s+(T 2) a /^ Xq(-^s) y 

2 Fg+Ca, VQg, s)e“^ s+<T 2) a / v 02 ^(-v^g^s]^ 

f 1 P 3± (-a, ( x,s)e< s ' w l) a /l J 


H dji 

Xq(-H,s)(h - v 02 ) 


+ 83.(8) P3±(-a,v 01 , S )e( s «l> a / v 0l 


V 


01 


x o( _v oi> s )( v 01 " v 02^« 


S € Sg-p (3*21) 


* {2 ks x^tt] 

F 2± ( a ^, S ) e - (s ‘ K, 2) a /l J x 0 ( S i,s)<I> lsv (4) 



+ S 2 (s)F £a; (a,v 02 ,s)e‘ (s+<T 2 )a/v 02 Xo(-v 02 ,s) y ^ 1 

* SW [/o' P 0d: ( - a '^ s)e(S+<ri)a/ " 


Xq(-H,s)(^ + v) 


+ 6 i( s ) F i ± ( * a , v Q1 , s ) e ( s+a l) a / v 01 ___ 


Xo(-v 01 ,s)(v 01 + v) 


0 < v < 1, 


(5-22) 








2k 


and 


J]+( v Ol) **%:(-a>Voi,s) | c 1 d 1 v 01 n , ls (v 01 )e"^ s ' Hy l^ a / v 01 


*(§ k s 


Xq(- v OP s ) 


5 


4 dp. 

H - V 0 l 


F 2±(^^ s ) e "^ S+CT2 ^ a ^ Xq(-M^) 

+ 5 2 ( s ) P 2 ± ( a > v 02> s ) e " ^ S+ff2 ^ a ^ V ° 2 Xo (- v 02 , s ) . V ° 2 — 3 

v 02 - v oi| 


Xq(-H>s)(h + v 0 i) 


+ | Pl± (-a,v 01 , S )e(B« 1 )a/v 01 


8 e S U . (5-25) 
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In the above equations, we have used the Xqj functions which Kuscer and 
Zweifel (ref* l4) have shown are continuous across the curves Gj in 
the s -plane (see App. A and Fig. 2) . For two material media, we take 
the ratio of their single -medium Xqj functions, 




(5.24) 


where 


Xoj(z,s) = 


(v 0J - z)X_) s (z), S 6 Sji 
(1 - z)Xj S (z), s e Sje 


(3-25) 


and 


tmmm m* 
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(5-26) 


For Re(z) < 0 , Xq(z,s) given "by Eq. (3*2^) is a nonvanishing analytic 


function of 
of v 0j (s), 


z 


and s provided s i £-o,j,- 0 j(l - Cj)J 
= 1,2. The quantity 


the branch cut 


k s * 8 ( 0 ^ - c 2 ag) + a x a 2 ( C 1 - c q) (5*27) 


is related to the difference between medium 1 and medium 2 continuum 
solutions; several equivalent expressions for k s are given in 
Appendix C. 

In Eqs. (3*15) we introduced the coefficients Ej + (v) since they 
are the forms of the normal-mode expansion coefficients which are 
extendable to the complex plane (refs. 2, 3) • Thus, Eqs. (3*1 6 ) through 
( 3 . 23 ) can be written in a compact form valid for Re(s) > -a m . These 
equations (see App. E) are 

%Xz,s) = l2±(z,s) 

± _ks_ q 2s(”) Xq(-z,s) r Eg i: (z’, S )Xo(-z', S )e' g(s ' H? g )a</zl ^ 
cgOg 2*1 i/(jt S!g s (z , )(z' + z) 


(5.28) 


EjJs.s) = I 1± (z,s) 

J- c 2 a 2 


fc s r 

* c 2 ct 2^<)( s) 2rti J fit 


E 2 ±(z * , s)Xq( -z ' , s)e‘ 2 ^ s+0 2) a / z 


«2s( z ’)( z ' “ z ) 


dz' , 


(3-29) 
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x ^ z ’ e) + {« 

x Lat(a,zSs)Xb(-BSs) J t 


c 2 a 2 fi 2s( z ')( z ’ + z ) 

22_g(°°) p L^+( "&,» z ’ ^ <i) dz 1 


l 00 / p 

^2s(°°) ^ c' M z ’ 


TTJ 


} 


(3.50) 


and 


Il±( z >s) = q; 


*0+( -a, z, s) e" 2 ( s+<y l) a / z ± ( a z s ) .£ Sa 2l 

c 2 a 2 S± v ' ' ' \0rfi]fc(-z,s)j 

CP L2f(a,z', 5 )Xq(-z»,s) 

\JJ c i c 2 ff 2^2sA z ' )( z ' - z ) 


%s(°°) P L 0+(- a > z S B ) dz’ ^ 

+ l 00 ) J c i ^lXol-zSs^is^'Hz' + z ) dZ J 


where for Ee(s) > -tjj, 

**<*■*•> -/* J) *' l * Wa)(Ma)/ “(ivi/ 0 


(3-31) 


0 f j±< x o’-^ jT77 


IL 

-2 c 3 a if 0 f 3±M rr 


(i - z 


+ ? f j±( x 0^ z )^js( z ^ ^O' 


(5.32) 


1(1) = -00 and 1(2) 


-a. 


(3.33) 
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In the above equations, z does not lie outside the contour C' which 
encircles as shown in Fig. 3 and % is defined as 

% h minCopOg). (5* 3*0 

The restriction Re(s) > -a m is discussed in the next section. The 
Lj+ functions were introduced as 

Lj ± (jc,v, s ) =P J± (x,v ) s)n+ s (v)nj s (v)e‘ (s ' w d )x/ ' 1 ', 0<V<1. (3.55) 

That Eqs. ( 5 . 28) -( 5 . 31) reduce to Eqs. (3.l6)-(3.23) as all 
contours C* are collapsed onto the branch cut v e ( 0 , 1 ) due to 
V z) can be seen as follows, if s e Sji , has a pole at 

z as. Vqj whose residue leads to a discrete term. When s e S^ e , 
fij S (z) does not vanish. The continuum terms are simply those due to the 
integration around the branch cut. 

The solutions fjc±( 3 C,p, s) and ^p+(x,(i, s) can now be written 
similarly as 



! 

S 
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2 rti 




f 

J n t 


E 3 +(z , / s)e^ s+<y l)( x+a )/ z * 
^ s (z')(z' + li) 


dz'j x < -a 


— r 

^ 2 jti iJ q i 


E3±(z%$)e^ s+a l^ 2 " a ^ zl 

/ • \ / » \ ^ j x > a, 

C* %s( z )( z - 1-0 


(507) 


for Re(s) > -a. 


m> 


t2p±( x ^^ s ) = 


{/, 


Lp|f (x^ z 1 , s) 
®W C « fi 2 s( z, )( z ' - n) 


dz' 


+ f La±(-x,z',s) J } 

"Jc’ n 2 s( z ')(z l + »j) J’ ^ 


< a, 


(3.38) 


for Re(s) > -erg* and 

r 

1 



Ia+(x,z",s) 

2 iti \ J P , «is( z, )(z' - H) 


dz' 


/ 

u n 1 


Mt(x,z',s) ± Ll^ ; (-a,z , ,s)e“^ s+crl ^ a “ x ^ z, 


%s( z ’)( z ’ + t*> 


dz'y. 


^lp±( x ^^ s ) = \ 


x < -a 


1 f T %(-a,^s)e‘( sw l)(^/ z + Mt(x,z',s) 

5tT vL, J2 ls (z , )(z’ - n) dZ ' 


1 / 
u P I 


L3±(-x,z»,s) 

C’ %s( z ')( z ' + n) 


dz * \ f x > a, 


(3.39) 


for Re(s) > -c^. 

The functions Mf(x,z,s) are also integrations over the initial 
distribution f l±(* ,|i) and are given hy 
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*<-■••> - ■/; wA ft •.*>// «**%»> ^ 


-k c i°i f f a±( x o^-^) 

«- \J r\ " ■* 


dn 

H - z 


+ | f a+( x o>- z ) ft Xs( z )j ^0/ x > a. 


^(X,Z )S ) =/ x ' a e-(B^i)(xo-x)/ Z 4 C101 ^ f3±(xo>(j) pSt 


- I c l°l/ f ]+( x 0’^) fTT 


z fl±(x 0 ,-z)fl ls ( z )j too. 


x < -a, 


(5-^1) 

for Re(s) > -cr^ and z not outside C'. Again, the discrete and 
continuum terms which appear in Eqs. (3.3), (3**0> (5*8) and (3*9) are 
due to the mvu\ and branch cuts of Qj s (z) which appear in the 
integrands of Eqs. (3*56) -(3* 59) • 


fa 


IV. PROPERTIES OF TRANSFORMED SOLUTION 


Analytic properties of as a function of s must be 

investigated "before we can recover the time -dependent solution 
ty(x,p,t) according to the inverse Laplace transformation given by 
Eq. (2.2). We need to know the behavior of \|r + in some right-half 
s -plane. Before looking at the details, * let us briefly review some 
results of earlier cited work in which Case’s method was used. 

In the previously mentioned work of Kuscer and Zweifel (ref. l4) 
and Erdmann (refs. 8, 9), expansion coefficients could be found ex- 
plicitly and this aided in the extraction of the s-dependence of their 
transformed solutions. They find that the branch cuts of v Q ^(s) are 
inherited by the transformed solution so that the integration contour 
of the inverse Laplace transformation must be deformed around these 
branch cuts. For the slab problem solved by Bowden (refs. 1, 4) 
expansion coefficients could not be found explicitly but theorems of 
Lehner and Wing (refs* 1 6 , 17) gave the analytic properties of the 
transformed solution in the s -plane. In that problem, the branch cut 
of Vq(s) is not inherited by the solution. Instead, the transformed 
solution has a finite number of poles at values of s, say Sq, 
which lie on the branch cut of v Q (s), that is, on the real s-axis. 
These poles contribute a sum of residues as the integration contour is 
moved to the left of them in the s -plane. Furthermore, in these 
previously solved time -dependent problems there is a real number, 
say 7 ^, such that the integration contour cannot be deformed into the 
region Re(s) < 7 -^ for arbitrary values of x. We expect the present 
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transformed solution to exhibit similar properties, that is, \|r + 
may not he analytic for Re(s) less than some number 7^ when x is . 

arbitrary while for Re(s) greater than 7 1 it should be analytic f 

except for poles and/or branch cuts. Such singularities probably 

occur where Vqj(s) ^ as i ^ ls ^ ranc ^ cut. | 

We first note that for arbitrary initial distributions f(x,|i), 

^ s+ (x,u) is not analytic for Re(s) < -cr m . This is true since each i 

of the inhomogenous terms I^ + of Eqs. (5.28) and (5.29) contains both 
L-, + and Lo+ as can be seen from Eqs. (5*50) and (5-51) and therefore, 

x- c. 

in general, is not analytic for Re(s) < -cr m , where ar m is given by j) 

8 

Eq. (5.3^). In particular, we note that for |x| > a, i|r^ + (x,ji, s) 
never appears to be analytic for Re(s) < - <7 . However, for special 
cases of material properties and initial distributions, ^2+( x ;^s) 
can be shown to be analytic for -cr 2 < Re(s) < -<j^ except perhaps for 
poles . 

We now look at the behavior of \J r g+ for Re(s) > -a m . Recall 

that the transform plane for the present problem must be taken as a 

superposition of two "single -medium" planes, that is, one for each 

material medium in the problem. The expressions (3*5)> (3**0, (5*8) 

and (3*9) for the transformed solution were not defined for seC . and 

outwardly appear to be discontinuous at seCh . However this is not 

the case. The complex representation of E^ + given by Eqs. (3.28) 

and (3*29) shows that such coefficients are continuous across the 

curves C . . Thus it is seen from the representation of \| f given in 

J sx 

Eqs. (3*56 ) “(3 • 59) that \Jr g+ is indeed continuous across the curves 
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It is convenient to introduce at this time the solution of the 
associated eigenvalue problem, that Is, the solution of Eg.. (2.13) 
subject to the boundary conditions ( 2.1k ) and (2.15) with f ^ + (x,p) 
Such solutions, denoted with a bar, have the form 


r — pi _ 

r= b l±^-v (^U/s)6i(s) +J ® 1 ±(-v)t ;i (^ v )( x ;^,s)dv, X < -£ 


\f v (x,n,s) + t (x,n,s) 8 2 ( 

L V 02 " v 02 J 


t s± (x,fi) = 


f Q B 2±( V ^2v( x ^' s ) ± ^2( | x | 


= ±t»l±\|f v (x,4,s)6 1 (s) ± J B 1± (-v)t lv (x,»Ji,s)dv, x > a. 


(4.1) 


where obviously Bj+ and b^ can be Obtained from the E^ + given 
by Eqs. (3«28) and (3* 2 9) in the case f ^ + (x,p) a 0. As we shall see 
later, the solution \|r + has poles at those values of s for which 
the associated eigenvalue problem has nontrivial solutions. In 
Appendix F, it is shown that as the slab thickness becomes very large 
this eigenvalue problem has only trivial solutions for Re(s) > -a p 
except perhaps on the branch cuts of ( s ) . When the slab thickness 

is not large, we still expect that if the eigenvalue problem, has 
nontrivial solutions for Re(s) > ~o 2 , they occur only when s is real. 
This has been proved rigorously using the method of Lehner and Wing 


(ref. 25) for several problems which can be obtained as special cases 

of the present problem: the bare slab considered by Lehner and Wing 

(refs. 16 , 17) and the slab surrounded by pure absorbers considered by 

Lehner (ref. 15) and Hintz (ref. 10 ). In all of these problems, there is 

no scattering in the reflector and, therefore, no branch cut of v^s). 

As already indicated, the X^(z, s) function inherits the branch cuts 

due to both v Q 1 (s) and v Q 2 (s) and these branch cuts lie on the real 

s-axis from -a. to -a (l-c ) and may or may not overlap depending 
<3 j d 

on the values of material properties. Note that c^ has been taken 
less than unity and this insures that the branch cut of v lies 
entirely to the left of s = 0 . In previously solved time -dependent 
problems, singularities of the transformed solution always occur where 
the Vq^(s) has branch cuts. Since the analysis of Appendix F indi- 
cates that for large values of the slab half-thickness, a, the singu- 
larities of . for Re(s) > -0 also occur where the vfs) have 
s± m Oj 

branch cuts, we will assume for all values of a that the singularities 
of \|r s± occur on the branch cuts of Vq^(s). In any case, we show 
that the only other singularities of \| r ., Re(s) > -0 which could 

S— in 

occur off the branch cuts of v^(s) are P°l es > vhose residue could 

readily be added to the time -dependent solution. 

In order to see the behavior of on the branch cuts of v qj( s ) 

we first look at if in the region seS A . For this region, 

si li 21 

the expansion coefficients are given by the equations (see App. G) 
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5 x..x ,, ks fi 2s(“) ( v 02 2 - d 2 ) hg (n) 
24 ' " 2 (v 01 2 - (X 2 ) g 2 <d) 


/ ha(V ° 2) ± |* g -'^ 8) + t B 2+ (v)h 2 (v) ^^--1 , 0 < (» < 1, 

\J* + v 02 M ‘ V 02 J 0 2 - v + ^ 


T °2°2 B p (0 l‘ a 2 )a / tl 

+ 5^7 B 2* (,l) 6 


Bl±(-d) + B 2* (|l) 

X X* 


+ ks fl 2s(°°) tl l(> a ) / h -g( v Qg) + h 2(“ v Qg) 

2 fi i s (°°) g 1 (n)V - v 02 H + v 02 


pi 20 (v) 

+ / (v)h 2 (v) — 

v 0 ** C 1 0 1 I 


T M- 0 1>*1± - h 2< v 02> * M- v 02) 


where 


+ ( v 02 2 “ v 01 2 ) J B 2±^ v ^ h 2^ 


X Ph (-cd) -(s+o 2 )a/« 


dv 

v 2 _ v 01 2 ' 


h p (o>) = cd -2fi- — _ e 
2 X 1s (-cd) 


M.) = 2i44 e (s+ °l )a/ “ 

1 o 2s W x 2s (-o>) 


«_,(»*) = d n$ s (d) flj s (d). 
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In addition the eigenvalue condition 


0 ■ \C ^ ( v)b 2 (v) 

V 01 * 02 V 01 v 02 J 0 


dv 

V + v 01 


( 4 . 6 ) 


must he satisfied. Since the eigenvalue condition (4.6) has different 
limiting values as s approaches the branch cut of Vq^(s), we conclude 
that there are only trivial solutions of the associated eigenvalue 
problem on the Vq^(s) cut. When s belongs to the branch cut of 
v Q2 (s), which is not also part of the v Q1 (s) cut, that is, when 
Re(vQ 2 ) = Im( V qi) = 0, it appears that nontrivial solutions of the 
associated eigenvalue problem may exist. From Bowden's results 
(refs. 1, 4) for the bare slab, it is expected that Eqs. (4.2) and 
(4.6) are satisfied only at isolated points, ( S n> . In the limit 
cgcrga -» “ these points lie on the branch cut of Vq^s), that is, the 
s n are real. The "thick-slab" eigenvalue condition is seen from 
Eqs. (4.2) and (4.6) to be Eq. (4.6) with B 2+ (n) = 0. 

If material properties are such that -a 2 < -a^, then a portion 
of the branch cut of v^ 2 (s) ' Lies in s€ ®2i^ S le* Iri ’ bhis re 6 ion 
however, s < -a m = and for such values, the solution \jr s+ (x, n), 

jxj > a, that is is not bounded as |xj -» °°. However, ijr may 

have nontrivial solutions on such a portion of the branch cut of 
V 0 2 (s). The equation for B^ and the additional constraint for this 
region are (see App. G) 
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B (n) a ± — (v 8 . m2) X 1 s(-H) 

21 2 E 1 S H 02 ‘ ei^IT x 01s (- M ) 


x J h g( v og) x ifl(“ v Qg) + h g(- v oa) x ig( v 02) 

V + v 02 X 01s^" v 02^ ^ “ v 02 X 01s^ v 02^ 


r 1 ^(v^v) Ais 


X la (-v) dv 


7 r — — — ) y 0 < j*t < 1 

(^v) V + (1/ ' - * ~ 


(M) 


and 


0 = h 2 (v 02 ) 


x ls(- v 02) + h / . x ls( v Qg) 

W-v 02 > ' 2 02 WW 



(it. 8) 


As we shall see later, the zeros of Eq.. (4.8) can, under some condi- 
tions, he poles of i|r 2± and therefore may contribute discrete modes 
in ty(x,n,t), |x| < a. For this reason we are interested in where 
these zeros lie and shall refer to them as pseudo-eigenvalues. 

We now show how the solution of the associated eigenvalue problem 
\jf + , is contained in the inhomogeneous solution, \|r s +, by following a 
procedure of Bowden and Williams (ref. 4). In Appendix H, it is shown 
that the original expansion coefficients of Eqs. (5*5) an d (5*4) can 
be written as 


y\\ 
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Aj±(fi) * jag± + A F 2 ±(a, v 02 ,s)J B^ ± (ji) + B J± (n) 


and 


a l± * + | %( a ;V 02 >s)J Bjj. + b 1± , scS^ A S 


2i (4.9) 


where B J+ and b^ + are given by Eqs . (4.2) - (4. if). The coefficients 


B^ + and b^ are given by 


B 2± (V) = ^S F l± ( - a ' V ' s)e 


(oi-a 2 )a,/v 


c 2 a 2 


ks fi 2sH (v 02 2 " y2) 


2 W (v 01 2 •• v 2 ) g 2 (v) 


(/o' 8 2t ()1 ) h 2^> 77 


. 1 w .. ^^2(^02) -r M" v 02)1 

+ ^ F^(a. v 0 o,s)|j n — + V *L v *j. j 


+ J" Q F 2±( a, ^ S ^ h 2^^ ^ + 


n + v 


+ X F a± ( - a '^ s)h i (,l) !'° 1 a ~ % ^^4 dtl> , 


v 02 * M c 2 a 2 v 






■■■■ 4 * 


39 


B (-v) + §2®£ B ;M .(v) e (<J:L ' C ' 2)a/V - ± 

1±' °1«1 3±' S g l(v) 


X< | F 2±( a »''o2' S> 



y h g( " V Q2 )~ 1 

v + v 02 J 


+ J Q {»2±<^) + F 2± (a, 4 ,s)jh 2 ( ( i) dd 

d 0 ( v q2 - f ) |i + v 

Ml 

[ CoCTo ( O'-] -CTp) a /' ; "l 

P 1 ±(-a,v,s) - P 2 ±(a,v,s)e J (4.11) 


and 


The coefficient 


^l^" v 01^1± * F +( " a v v oi' s “ Pl±* 
|*2± + 2 F 2±^ a,V 02' s ^j is 6iven hy 

r . i „ / n ’ v oiPi± + Pa± 

L J ( v 01 a l± ” a 2 ±^ 


In these equations, <?^ + and 0^ are 


(4.12) 


(4.13) 


a l± = h 2^ v 02) ± h 2^‘ v 02) + ^ v 02 2 " v 01 2 ) f Q B 2i(M)h 2 (a) 2 , 

ol_ 

a 2± = v 0sM v 0g> T v 02 h 2^ v 02) + < v 02 2 - V 0 1 2 W B 2± (n)h 2 (a) —mL. 

0 * ~ V 01 


(4.14) 


4o 


Pl± - | F 2±< a ' v 02’ s >[ h 2< v 02 ) + h 2 ( - v 02 ) ] 


+ (v 02 2 - v 01 2) 


+ P 2± (a,(i,s)jh 2 (n) -g-2i 


± (v 2 - Vr,r} 
V v Q1 Vq 2 


*)/ o F^-a^s)^^ 


2 2 
H “ v 0 2 


± [f 1+ ( -a, 


V 01 , S )h x ( v Q1 ) + F-jj, ( -a, -v 01 , s )h 




&2± ~ 2 F 2±( a ' V 02 f8 ) v Q2 


Jh 2 (vo 2 ) ± 


h 2 (-v 02 ) 


+ ^ V 02 2 " v 01 2 > / B 2+(^) + %fc(a,li,s) h 2 (n) 


2 2 
* - V 01 


+ (v, 


2 - v 2 
L v 02 


p 1 

) / F 1+ (-a,fi,s)h 1 (n) 

u r\ 


V? “ V, 


Ol'li:' “ a> v 01> ■ n i' v Ol/ ’ Wl±^^- V 01^^ n l^- V 01' 


(4.15b) 


kl 


In terms of these quantities, the solutions can he written as 

t 2± (x,|i,s) = ja^j. + | F^Ca^Vo^sJj^C^^ 8 ) 


f Q B 2fc( v )^2v( x ^> s ) ± *2(-v)( x ^> 6 )] 


dv 


+ j jF 2 ± (x,v,s)t 2 v (x,M,s) ± P 2 +(-x,v,s)'|r 2 (_ v )(x,n,s)Ja,v 


+ |jr 2i (x,V 02 ,s) ± IVyX -X, -V 02 , S )Jy v ^(x )( l, S ) 


|[f 2± (x,-> 


v 02> s ^ “ ~ x > v Q2 fS 


) U.v ( x #^s), | x |< a, 
J v 02 


(4.1 6) 


and 


* 


2+(x,n,s) = fa2± + ^ F2+(a,Vo2> s )j ♦l±C 3t >^ B ) 


±|hi+ - P ± (-a,v 01 ,s) + F 1 ± (-x,-v 01 ,s)J\|f Voi (x,^s) 


± F 1 ± (-x,v 01 ,s)\|r_ v _ (x,^,s) 


01 


± r[- ( - 


v) - F+(*a,v,s) + Fi+( 


-x, -v,s)J^ lv ( 


x,|i,s) 


± / Pa±(-^v, 8)l|r l{-v)^^ 8 > 4v ' x > a ‘ 


(4.17) 










—.u. 




The solution tx±( x >H> s ) for x < -a has a similar form. In these 
equations, \J^ + (x,n,s) are the parts of i|f g+ (x,n) which are given 
by Eq. (4.1). Equation (4.4) is written in terms of a^j. as 


T h 1 (-v 01 )b 1± - o^. 


(4.18) 


Consider now what happens on the branch cut of v Q ^(s ) where 
Vqi * i | v oi| f° r * m ( s ) = 0" and Vq 2 = -i |v Q ^| for Im(s) = 0 + . 
From the above equations it can be seen that the quantities B^, B-jj., 
B^, Bj 4 _, Oj+, ag+, Pxt an< i $2± n °b inherit the branch cut of 
Vqi(s), Equations ( 4 . 18 ) and (4.12) show that b^ and b have 
branch cuts due to that of Vq 1 (s). Equation (4.15) indicates that 
[ &2t ♦ | F 2± (a,v 02 ,s)J has the branch cut due to Vq-^s) unless 
<3t, l±/ a 2+ is e< l ua b bo 3 i+/ 32±* g enera b; this will not be true since 
Pl+/P 2 ± depends on the arbitrary initial distribution f+(x,n) whereas 
C! l+/ 0( '2t does n °b* Therefore, it is concluded that both and t p+ 

inherit the branch cut of Vq 1 (s). 

On the branch cut of Vq 2 (s)> the quantities B^, B^*, b^., 3-^+ 
and 3 2 + are single -valued. Since the quantities and of 

Eq. (4.l4) are related above and below the branch cut of Vq 2 (s) by 



it follows then from Eq. (4.15) that on that part of the branch cut 
of Vq 2 (s) which is not also part of tv, Vq^(s) cut; that is, for 
ReCvQg) * Im(v’oi) = 0; we have 






(4.20) 


jag± + | F 2± (a,v 02 ,sjj + = ifa a + | F 2± (a, v Q2 ,sjj 

if the denominator on the RHS of Eq. (4.15) does not vanish. It is 
seen from Eqs. (4.1), (4.2), (4.5), and (4.4) that for this same region. 


[t s± ( x ,ii)J = ± [, S ±( X ,)] . 


(4.21) 


Hence the product 


j^2± + i F 2±( a ' v 02’ s )] 


(4.22) 


which appears in t s ± does not inherit the branch cut of Vq 2 (s). 
However, the denominator of ja^ + 1 F 2± (a,v Q2 ,s)J , namely 

( V 01 a i± " a 2± ^ :i ' s e< l uivalen ' fc to the eigenvalue condition, Eq. (4.6) , 
Thus, if the associated eigenvalue problem has a nontrivial solution 
at s = s n , Re(s) > -o^, then ^ s+ has a pole there. 

We briefly summarize the analytic properties of the transformed 
solution t g+ (x, n). For arbitrary initial distributions f + (x, fi ) , 
is not analytic to the left of Re(s) = -a m in the s -plane, 
whereas to the right of Re(s) = -cr m it is analytic except for the 
branch cut along f (due to the branch cut of v^(s) ) 

if a m > a^(l-c 1 ) and poles at the values of s at which the 
associated eigenvalue problem has nontrivial solutions, \|r + . We have 
assumed that for arbitrary slab thicknesses, a, these poles, if they 
exist, lie on the branch cut of Vq 2 (s) since this is the rigorous 
result obtained by others for several special cases of the present 


problem and. obtained herein for the case when CgCfga is large. For 
special values of material properties and initial data, f , (x.fi) for 
| xj< a (that is tg+) way *> e analytic in the region -erg < Re(s) < -a^ 
except perhaps for poles. 











V. RECOVERY OF TIME -DEPENDENT SOLUTION 


The time -dependent solution t( x >H>t) is obtained from the 

inverse Laplace transformation Eq. (2.2), where 7 is to the right of 

all singularities of t s (x,n) in the s -plane. From the analysis of the 

preceding section we expect that we can choose any 7 > max ^l(^*“®2.) > 

-0 2 (1-C 2 )j . In order to show the time dependence of the solution 

ty(x,n,t) more explicitly, we deform the inversion contour as far as 

possible to the left in the s -plane by making use of the analytic 

properties of \|r (x,n) obtained in section IV, 

We first look at the behavior of i|r s+ on the contour Re(s) = 7 • 

This contour crosses both of the curves C. and it has been shown that 

J 

is continuous across these curves. As | s| -* <*> on such a contour, 
seSie A S 2e and we show in Appendix % that \|r + behaves as follows 


ox -(s+ 0 2 )(x-x o )/n r n r\ 

t 2 ±(x,H,s) -+-J e f 2±( x 0'^ + °(l)J dx O 


-(0 2 -0l) 


'(n) P-a -(s+0 1 )(^^ 


jjl±(x 0 ,M) + 


(5.D 


for |x| < a and |i > 0; 


tof( x >^ s ) 




^>(w) 


■ (s+al) (^ 


f i+( x o>-M) + 


(5.2) 
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for |x| < a and p < 0; 






/ \ 2a 


>a -(s+or 2 ) 




f 2±( x 0'^ + °(l) K> 


-a -(s+a.) 


1 ^ / |f 1± (x 0 ,n) + 0(1) dx 0 


for x > a and p > 0 and 


r»X -(s+a^f^p)*- /n\“ 

+ e v [ f i± ( v^ + °(i)J to o 


t 1± (x,^s) jlj -J e ^ S+<? ( |nlj f^CxQj-lxl) + o(~j dx 0 (5^ 

for x > a and p < 0. Expressions similar to Eqs. (5*5) and (5*^-) 
are obtained for x < -a. It is seen that t s + is not necessarily 
0^1^. However, the parts which are not can be easily inverted as 
follows. Define for all s the function f^+Cxjii) as that part of 
each of Eqs. (5*1) - (5**0 which is not °(|). We show in Appendix I 
that upon making the substitution 


x - x Q = »it, p > 0 

x - x = |p|t, p < 0, 


that \lr (x,p) can be written as 
us± 


♦us ± (x, '* ) e [w*^] dt ' 


That is, the parts of t s+ vhich do not "behave as as j s J 

Re(s) * 7 can be inverted by inspection. The solution 
is given by 


r -pgt 
e f^Cx-fat^) 


, t < 


a+x 


^(x^t) = ( 


" a lt **( cr 2 " cr l^ *£* a+x 


e e 


(5.7) 


for | x | < a and |i > 0 ; 


'l ; u +( x ; M^t ) 


" -o 2 t 

e f , 2±( x "^ t ^) 


t < f “ x 


M 


(02-01 ) — jr~ 

e e 11 f 1 ± (x-nt,n), t > £=£, 


for | x | < a and U < 0 ; 


= < 


■it 

e f 1± (x-Mt^) 

-a 2 t (d 2 -ai)^j~^ 


, t <££ 

9 n 


(5.8) 


f2+(x-|it,|i). 


x-a 


— <t < 2 ±a 
H M 


( 5 . 9 ) 


-CJnt -(Up-^) 

4 f 1 ± (x- 4 t,n), t > £±a, 


e e 

v- 


for x > a and 4 > 0 and 


-a-jt 


♦u+(x>M,t) = e fj+tx-nt,^) 


( 5 - 10 ) 


for x > a and n < 0 . That *u± describes the motion of uncollided 
neutrons from the initial distribution can be seen by direct 
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substitution, that is, satisfies the uncollided equation 

d\lr , , 

M + tf(x)t«fc = 0. 

In the limit t -> 0, we note that 


= f ± (x,|l). 


( 5 - 11 ) 


For arbitrary f(x,|i) which vanishes as 
given by Eqs. (5.6) and (5-7) - (5*10) is an analytic function of s 
for Re(s) > -cr m for almost all x and U. If f^ + ■ 0 (fry. = 0), 
then y us+ is an analytic function of s for Re(s) > -a 2 (Re(n) > -cri) 
Therefore the function <J> s+ (x,u) defined as 

® s± (x>n) ■ ^ s± (x,u) - t us± (x,|Ji), Re(s) > -o m , (5.12) 


has the same analytic properties as i|r in the right-half plane 
Re(s) > -o m except that it is as |sj-*«>. If has a 

branch cut along jja m , -^(l-c^ ; i.e,, if a m > or 1 (l-c 1 ); then 


[*-] - [*s ± ] + = v] - *s± + - 


(5-15) 


Similarly if t s ± has a pole at s = s n , then 


Residue 


M Sn = Residue Ws n - 


( 5 . 1 ^) 


The definite parity parts of the time -dependent solution therefore 


can be written from Eg,. (2.2) as 


A/+1W , 

t+(x,n,t) = + gjg J <I> s± (x,n)e s ^ ds. (5-15) 

T i* 

Now using the analytic properties, we can deform the contour to the 
left and obtain in general 



+ Si p^O ]„ ' tr s±< x ^> eS ' fc iB ’ -°m < -Old-' 1 !) < v 

C P 

(5.16) 

where Cp is a small circular contour of radius p with center at 
s = -cr^(l-c^) . Generally the point s = -<j^(l-cj.) not satisfy 

the eigenvalue conditions, Eq. (4.6), and the contribution from the 
contour Cp vanishes as p -» 0. If however s = -a^(l-c^) happens 
to satisfy Eq. (4.6), the contribution from the contour C p has the 
form of a discrete residue term. Details concerning this point are 
discussed in Appendices I and K. 


Equation (5.16) is the solution of the time-dependent problem 
written in a form in which the uncollided portion of the initial 
distribution f(x,n) has been separated, For arbitrary f(x,n) the 
contour cannot be deformed further to the left. In the final section 
it will be shown how this solution reduces to those obtained previously 
by others for special cases of the present problem. 

We close this section by indicating the form of some parts of 
Eq. (5.16). The uncollided term, ^(x^n^t) is given explicitly by 
Eqs . (5.7) - (5.10). The form of \|r s± (x,n) the branch cut 

“ a l( 1_c i)J was 8 iven in section IV. From those results, it is 
seen that on this branch cut [.„< X,n)l - j^ s± (x,|i)J can be written 
from Eqs. (4.16) and (4.17) as 


^ - [*s±( x '^ = ^[ a 2± 


|^ a 2± + I P 2t^ a ' v 02 ,S |J 


- |a,2+ + — P p±( a , v 02^ s ) r) ?2±( x '^ s ) 


for |x| < a and 


[* s± (x„j] - [ V ^j| = Ja^t + | P2t(a,v 02 ,s)j u, s )J 

•|®at + I F 2±( a ' v oa’“|J |^i±( x »w,s)J 

1 |*i* - F ± ( - a ' V' s >] v (x ' M ’ s ’ 

- fb-i+ - F+(-a, V(y.,s)j ♦ _ v ( x jH>s) 


(5.18) 
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for x > a, where a^j. + | 1? S!:( a ^ v 02' s ) is Sl ye ^ by Eq,. (4,12), 

. , V». ^ #v n 

^ 4 .(x,u,s) by Eg. (if. 1) and b 1± - F+( -a, v Q1 , s ) by Eg,, (if. 12), The 
solution V 5+ (x,|i) bas poles at s » Sq, . . .,s^ due to the poles of 
a 2± + | F 2±( a ' v 02' 8 * VQgC^l)/^^ Again, from the results given in 
section IV, it follows that 


Residue 


|t s ±(x,|i)e B *j 


*=+(*. 


,n)[v 


-( m )/2 


X Residue <v 


(1TXJ/2 


J®2fc + \ F 2±( a > v 02> s ) ^ 


Note that the factor v. 


(it 3) /a 


is introduced so that 


and a# + | ^(a, v Q2 , s ) v 02 


(lTl)/2 


^s± v 02 


(5-19) 

-(l+l)/2 


are single -valued on the branch 


cut of vq 2 [cf., Eels. (if. 20) ; and (4.2l£) . These terns have an 


exponential time dependence, e , and we have obtained the implicit 


Eqs., namely, (4.2) and (4.6), from which the eigenvalues On) can 


be computed. Since information concerning the behavior of eigenvalues 


(i.e., number, location, etc.) as a function of material properties is 


not readily obtained analytically from such expressions, we have made 


a numerical study of real time eigenvalues and the results are discussed 


in the next section. 


-1 




2 


VI. CALCULATION OF BEAL TIME EIGENVALUES 


We firs>5 note that the eigenvalues and pseudo-eigenvalues depend 
on five parameters (c,^, c 2 , dg> and a) and therefore many numerical 

computations would be required in order to see the specific dependence 
on each parameter. As we shall see, the bare-slab results of Bowden 
(refs* 1 , 4 ), the theorems of Hintz (ref. 10 ) for slabs surrounded by 
pure absorbers and Some observations of the present numerical results 
for a few reflected slab cases allow us to draw some conclusions about 
the behavior of eigenvalues for reflected slabs as a function of the 
slab half -thickness a. However, rather than compute eigenvalues ^s n ^ 
in terms of cq_, dp, Cg, dg> and a we define a nondimensional 
variable £ and nondimensional parameters d R , dp, and A as 

i = ** “ W 0l) = cgOg and A = °2 a 2 a - (6>1) 


In terms of these quantities, the branch cut of Vpg becomes the real 
interval ( 0,1) and the branch cut of v Q p becomes the real interval 
(-dp, -dp + dp) • Since d^ and c^ are non-negative, it follows that 




( 6 . 2 ) 


where the equality holds only if dp = 0. Also we have restricted 
Cp < 1 so that -dp + dp ^ 1 implies that Cg < 1. Obviously, 
dp = 0 when the reflector is a pure absorber or a vacuum and dp = 0 
when the total macroscopic cross sections of the two media are the 
same. We have seen from the last section that in general the inversion 
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contour can "be deformed to the left only as far as Re(s) = -cr m , which 
corresponds to Re(£) = max(-a D ,0) . However, there are no eigenvalues 
on the branch cut of v Q1 so the region of the real £ axis where the 
eigenvalues &> should appear is 

max(-<jp + < tn < 1» (6.3) 

This interval corresponds to s e S 2 £ and Eqs. (4.2) and (4.6), 

written in terms of the quantities of Eq. (6.1), are solved numerically 
to obtain the real eigenvalues for specified a R , and A. 

In addition the pseudo-eigenvalues are obtained numerically by solving 
Eqs. (4.7) and (4.8) also written in terms of the quantities of 
Eq. (6.1). Numerical results are also obtained in the thick-slab 
approximation, that is, Eq. (4.6) with 5g+(|i) = 0. Details concerning 
numerical procedures and computational equations are given in Appendix J. 

The calculations were done on a Control Data 6600 computer system 
at NASA Langley- Research Center. 

The time dependence of discrete modes is seen from Eqs. (5.1 6) 
and (5.I9) to be 

e Snt =? (6.4) 

Now £ n m -(jp + a R implies that s n * -<ji(l - ci) < 0 since c-^ < 1 

and the equality holds only if cr^ = 0. Therefore such ^ correspond 

to time-decaying modes regardless of the value of c 2 . For values of 
^ within the interval (6. 3)> the time decay or growth depends on 
whether c 2 £ n is less than or greater than unity as can be seen from 

Eq. (6.4). A discrete mode represents a critical system if c 2 £ n = 1. 


■s* 
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The largest eigenvalue £q with an even parity eigenfunction corresponds 
to a critical slab problem with 


and 


c 


1 

slab " 


c reflector 


gR 

£o + ff D 


^slab a critical “ 


( 6 . 5 ) 


where ^critical ^ ne critical slab half -thickness. For a bare 
sphere (a^ = 0) the largest eigenvalue ^ with an odd parity eigen- 
function gives the critical sphere radius, a C riticai> vhen i-b is use<i 
in Eqs. (6.5) in place (ref. 21). 

Many different combinations of material parameters could be 

f 

considered, but here we restrict our study of the eigenvalue behavior 
to the case of overlapping branch cuts. As 0 ^ departs from zero, we 
would like to see how the eigenvalues depart from those previously 
reported (refs. 1, 4) for a bare slab. A comparison of the present 
eigenvalues for vacuum reflectors, i.e., a# = 0, with those of 

Bowden (ref. 1) is given in Tables I and II. Results generally agree 
to three figures for slab half -thicknesses A from 0.4 to 20. In 
Table II, eigenvalues calculated in the thick-slab approximation are 
also shown for bare slabs. For slabs with half -thicknesses A > 1, the 
thick-slab approximation generally agrees with the numerical solution 
of the exact eigenvalue condition to three figures. This can be seen 
from Table III where we compare such results as Oj| departs from zero 
with ffp - 0. From the bare slab results (ff^ = 0) of Tables I, II, 
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TABLE II.- EIGENVALUE 


£ 0 FOR THIN BARE SLABS 


Slab 

thickness 

A 

Thick-slah 
approximation 
Eg.. (J. 19) 

Present 
Eq. (O’. 6) 

Bowden 
(ref. 1) 

1.0 

0.702 

O.705 

O.705 

.8 

.612 

.615 

.615 

.6 

A75 

.485 

.485 

.4 

.241* 

,282 

.282 

.2 


.045 

00 

O 

• 


*No solution found for £ > 0.001 





III. - EIGENVALUES £n MD t, FOR THIN REFLECTED SLABS 
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and III, critical slab half-thicknesses are obtained from using 
Eqs. ( 6 . 5 )* These are compared with critical slab half -thickness 
results of Mitsis (taken from ref, 7) in Figure 4 (open symbols). 

Closed symbols give critical sphere quarter-diameters obtained from 
Eqs. (6.5) and £1 while Mitsis' critical sphere results are taken 
from ref. 20. The agreement is good to the scale of the Figure, For 
Cj> = 0 the eigenvalues and ^ have also been compared directly 

with numerical bounds computed by Mullikin (ref. 21) for bare slabs 
and spheres and again the agreement is good. Critical half -thicknesses 
of slabs with infinite reflectors have been recently computed by 
Kowalska (ref. 12) for a number of combinations of Cg^t, and 
c ref lector* Some Present results £q for 4 0 can be compared 
with her critical slab half -thicknesses. Her parameters are given 
in terms of and present input quantities a-Q, and A by 

Eqs. (6.5)* Figure 5 gives a few present cases (circles) for which 
c slab was close to some of Kowalska's points (diamonds) (ref. 12). 

No attempt has been made yet to compute points which lie on Kowalska’s 
curves. The present cases for c p -| ^ ~ 1.11 are from A = 2 and 
1.4 in Table III. 

The remainder of the results have been computed for A = 5* For a 
bare slab with A = it can be seen from Table I that there are five 
eigenvalues. We have studied the behavior of these eigenvalues as 
Op departs from zero for several values of < 7 p. In Figure 6, results 
are given for <j d = 0. Our calculations show that the largest eigen- 
value, is present up to <x R = 0.9999* Apparently this eigenvalue 

remains up to crp = 1, which is only obtained for c 2 < 1. All other 


critical dimension 
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Critical slab half-thickness (ref. 7) 
Critical sphere quarter-diameter (ref. 20) 
Critical slab half-thickness Present 

Critical sphere quarter-dlamster (r 

O and • to scale of figure 


Mltsis 


Present and Bowden, 
(ref. 1) 


Figure 4 .- Critical dimension of bare systems. 


•• 







62 


eigenvalues disappear into the branch cut of Vq^_ at £ n * cr R , labeled 
with *, which corresponds to a time-decaying mode, regardless of the 
value of Cr,. In fact, on Figures 6 - 9> we indicate the points at which 
an eigenvalue or pseudo-eigenvalue coincides with the branch points of 
Vq^ by an asterisk, *. Even though such points appear to have a 
discrete eigenvalue type of time dependence, we feel that they are 
properly part of the branch-cut integral contribution. We note that the 
branch points of are located at £ = -a D and £ = " (T D +0 R and - find 

that the limiting form of the condition which determines whether or not 
such points are eigenvalues (or pseudo-eigenvalues) no longer depends 
explicitly on Or or a R . (See Appendices J and K. ) The theorems of 
Lehner (ref. 15 ) apply for o R = 0 in this Figure. 

In Figure 7/ results are presented for <jp * -0.65 + 0.5 <J R . These 
represent what happens for -Op in the range between zero and 


K n l } where the notation 
L 0J O R =O 


cr R =0 


means bare -slab eigenvalue, 


which we note depends on Cg,^ ®nd a * 1116 °P en and- closed circles 
represent eigenvalues as in Figure 6 while the half -closed circles are 
pseudo-eigftnvalues corresponding to s < -cr m = -oi» Again the largest 
eigenvalue, £q, appears to remain provided that Cg > 1* Here, as in 
the next two figures, results for o R = 0 agree with the theorems of 
Hints (ref. 10) which apply only for c-j_ » 0. Basically his result is 
that the strip Re(£) between 0 =lnd belongs to the continuous 

spectrum and that the bare-slab eigenvalues lying in this interval are 
not eigenvalues of the slab surrounded by perfect absorbers. He finds 
that there are no eigenvalues if -cfp > but says nothing about 
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the physical significance. It is seen from Eq. (6.2) that for such 

cases < l/c 2 and corresponds therefore to a time-decaying mode. 

In other words, stationary (critical) or time-increasing modes cannot 

disappear into the continuous spectrum as material properties are 

varied. In fact, we have seen that when <y R / 0 such modes cannot 

disappear into the branch cut of either. In Figure 8, results 

are given for -Op + a R = 1 which we recall implies Cg < 1» For this 

case, all of the bare-slab eigenvalues lie in the continuous spectrum 

found by Hintz (ref. 10) when cr R = 0. In both Figures 7 and 8, 

s = -(j m corresponds to £ = -dp. Figure 9 shows the behavior of the 
eigenvalues for = 1 and it is seen to be similar to that of 

Figure 6. For ff R = 0, the continuous spectrum found by Hintz (ref. 10) 

lies in the strip -a D = -1 < Re(£) < 0. Here s = -a m corresponds to 

5 = 0. 

All numerical results indicate that- real time eigenvalues &> for 
material reflectors are finite in number and tend to eigenvalues pre- 
viously obtained for a vacuum as or R -» 0, as do the pseudo-eigenvalues 
for s < -a m . If the set is empty, then the neutron density is 

necessarily decaying in time. Conversely, if the neutron density is 
stationary or increasing in time, then the set is not empty. 

One also expects that if Cg > 1, then a critical thickness can be 
found. That is, the largest eigenvalue £q must be present for large 
enough slab thicknesses for the given Cg. This can be seen from 
Table I as follows: For example, if -<y R = 0.8, then the eigenvalue 

for A = 1 is not present, while that for A = 5 would be and 
represents a mode whose amplitude increases exponentially with time 








and branch- 
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for c 2 > l/0.9T5. That for A » 20 needs only c 2 > l/o.998 in 
order to represent a time -increasing mode. 

As pointed out at the beginning of this section, we can draw some 
conclusions concerning the behavior of ^ or reflected slabs as a 



function of the slab half -thickness a. That is, given Cj and Cj 
what can be said about {in} as a function of a. We base the 
following conclusions on the observation that if at <Tr = 0 lies to 

the right of -a^, then it appears to remain to the right of -crp + Op 

as dp increases until -dp + dp = 1. (See Pigs. 6 and 7*) The 

dependence of on slab half -thickness is given in Tables I-III 

and many more points are given in reference 1. First, if -ff_ + d 0 > 1 

(recall that this implies eg < 1) then the set is empty for all 

a. However, there may be pseudo-eigenvalues if -dp > 0. Next, if 
-ffp + dp < 1 then two cases arise, depending on the value of dp. 

(a) When -dp > 0 then regardless of the value of Cg, we can find an 

a* such that a < a* implies that the set &> is empty, whereas 
a > a* implies that the set is not empty. The number a* is 


obtained from the bare-slab result 


0^=0 


|fo( c 2>02> a *)J 0RSsO * -Op* (6.6) 

(b) When -dp < 0, the set {in} is never empty, Thus, given c^, oj, 
a and the bare-slab eigenvalues corresponding to eg, dg and a, we 
can say whether or not the set is empty. Furthermore, the number 

of eigenvalues will not exceed the number of bare-slab eigenvalues 

fcu-o} which are greater than -a^. Finally, the number of real 
reflected-slab eigenvalues and pseudo-eigenvalues does not exceed the 


the number of real 


number of bare-slab eigenvalues. 


aSPJBKt! 



VII. CONCLUDING REMARKS 


It has "been shown using Case's method that the solution to the 
initial-value problem of monoenergetic neutrons migrating in a finite 
slab (properties Cg, arf with infinite reflectors (properties c 1; 
a-j) can be written in the form 

t(x,[i,t) = \Jr u (x,n,t) Residue 


+ C1) ([ + S (x ^ ) ] - (> ( H + 


st , 
e ds 


•'-ff m +ioor -i 

|'i r s( x ^M) “ t us (x,n) eS ~d m < " c l) ^ s n* 

-CT^-ioo L J 


(7.1) 

In this equation, t is the neutron speed multiplied by the real time, 
the minimum of and Og and each ijr function is the sum of 

its definite parity parts ijr^. Some terms of the solution (7.1) will not 
be present if -a m i -a x (l - erf i s n . That is, if -cr^l - Cl ) < -a m 
then the branch-cut integral does not appear. Likewise, if all 
s n < -ffj(l - C]_) , then there are no residue terms. These discrete 
eigenvalue terns are characteristic of a finite slab (refs, 1, k) while 
the branch-cut integral term is typical of a serai-infinite medium 
(ref. 14). The term i|r u (x,jx,t) describes the behavior of neutrons from 
the initial distribution, f(x,|i), which have not suffered a scattering 
collision and its definite parity parts are given in Eqs. (5.7) -(5.10) . 
The discrete eigenvalue terms in Eq. (7.1) are given by Eq, (5. 19) while 

the integrand of the branch-cut integral is given by Eqs. (5.17) -(5.18) . 
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The definite parity parts of the last integrand are given by Eqs. (5*6) 
and Eqs, (3»56)-(3»39) • The eigenvalues ( 3 n> can be computed as was 
demonstrated in the last section; thus, everything which appears in 
Eq. (7-1) can be readily calculated. 

In all special cases of the present problem which have been solved 
using the Lehner-Wing technique (refs. 10, 15 l6, 17), c-^ = 0. In 
these cases, there Is no branch cut due to Vq^(s) ; therefore the 
branch-cut Integral is not present in Eq. (7»1)» It was shown that as 
c-|_ -» 0 the eigenvalues, which are greater than -<r m approach 

those for a bare slab as do the pseudo-eigenvalues for s < -<r m . The 
solution ^ 2 + has the proper behavior as c^ -* 0 since those terms of 
Eq. (3*30) and (3*31) which appear to blow up in such a limit actually 
cancel when the contour C' is collapsed onto the portion of the branch 
cut of a Js (z'), 0 < z‘ <1. When the uncollided term is combined 
with the last integral it is then seen that the solution (7*1) and the 
eigenvalues {s n ^ have the behavior required by the theorems of 
Lehner (ref. 15) and Hintz (ref. 10) . The present problem reduces to 
those considered by Lehner and Hintz when 

c-^ = 0, o-j . ~ a 2> k 6 * 11161 * (^ef. 15) 

and 

c^ =0, ^ Oq; Hintz (ref. 10). (7,2) 

Hintz shows that for cr^ - his spectral results reduce to those of 
Lehner, 

In order to describe the same physical problem in the slab as that 
solved by Lehner and Wing (refs, 1 6 , 17) we must not only have 
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but also 


* 0 and cr^ = 0, 


f(x,n) * 0, 


rx < -a, 
Vx > a, 


H > 0 

\x < 0. 


( 7 - 5 ) 


( 7 - 4 ) 


In other words, neutrons from the initial distribution outside the slab 
cannot impinge on the slab faces at times t > 0. Lehner and Wing 
solved the time -dependent problem with boundary conditions 

t(±a,n,t) = 0; (1$0, t > 0. (7.5) 

Restrictions (7*5) and. (7*4) in the present solution make I^(u,s) 
and therefore s) depend only on slab properties. Then, in 

looking for solutions inside the slab (|x| < a), the inversion contour 
along Re(s) = -a m can be deformed back to Re(s) * -o 2 , and we pick 
up a residue contribution from any pseudo-eigenvalue in the region and 
thus obtain the Lehner-Wing results. That is, the solution has the 
proper form and all bare-slab eigenvalues are recovered. Hintz 
(ref . 10) did not indicate how the Lehner-Wing solution for the bare 
slab could be obtained from his results. Here we emphasize that he is 
not solving the same physical problem inside the slab unless both 
conditions (7*3) and (7.4) are satisfied. 

The analogous problem for c^ 4 0 in which the inversion contour 
can be deformed to the left of Re(s) = -<J m for |x| < a is obtained 
when <j 2 > and f-^(x,p) = 0. That is, if 

f(x,u) = 0, I x I > a and Ug > 0]_, (7*6) 


o 
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then all terms in Ig+(n,s) vhich contain (s + c^) factors in the expo- 
nentials are identically zero and this allows us to deform the contour 
along Re(s) = -a ra hack to Re(s) = -oq when |x| < a. Such a 
deformation is not possible for |x| > a; for this latter range of x we 
must stop at Re(s) = -0 m = -c x . If there are pseudo -eigenvalues in 
-Cg < Re(s) < - 0 ^ = - 0 m (see, for example. Fig. 7) they will appear in 
the solution for |x| < a as residue terms which have the exponential 
time dependence. They are not eigenvalues for the reflected slab though, 
since such terras do not appear for |x| > a. Erdmann (refs. 8, 9) 
solved the time -dependent problem for two semi -infinite media where an 
isotropic pulse of neutrons was introduced at the interface and found 
that the inversion contour for x e medium j could be deformed to the 
left as far as Re(s) » - 0 j . In the present problem, such deformations 
can be made only when conditions (7*6) are satisfied. It appears that 
the contour Re(s) * - 0 m cannot be deformed to the left of 
Re(s) » - 02 > since the implicit equation which determines A 2 ±(p, s) 

(see Eq. ( I. 3) ) requires Re(s) > - 02 * Apparently Re(s) = -0 2 is the 
edge of a continuous spectrum in all cases for the reflected slab. 

We briefly summarize the results which have been obtained. The 
present solution has been shown to have the required properties in all 
special cases which have been solved previously by others using the 
Lehner-Wing technique. However, in all of these rigorous solutions, 
there was no scattering outside the slab. We have seen that with infinite 
reflectors on the slab and neutrons anywhere outside the slab initially 
that it is possible for some neutrons which have spent their entire 
history in the reflector to impinge on the slab faces at times t >0. 
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Such neutrons have a collision rate which is characteristic of reflector 
properties and this, in general, restricts us from deforming the 
inversion contours to the left of Re(s) » -a m . We have illustrated 
two cases in which a further deformation is possible for |x| < a, by 
eliminating neutrons outside the slab initially which can later impinge 
on the slab faces. This is equivalent to a further restriction on the 
Hilbert space which has been used in some of the above-mentioned 
rigorous solutions. The exact eigenvalue condition has been obtained 
and real time eigenvalues have been calculated for a number of 
combinations of material parameters. The largest eigenvalues have been 
shown to agree with criticality results of others. Our calculations 
show that eigenvalues can disappear into the branch cut or continuum as 
material properties are varied and we point out that all such 
disappearing eigenvalues correspond to exponentially time-decaying modes 
regardless of the value of c 2 since we have taken c 1 < 1. We expect 
(but have not shown) that there is no drastic change in the shape of the 
solution given by Eq, (7.1) when this happens; we conjecture that one of 
the integrals in Eq. (7»l) probably has resonance -like terms due perhaps 
to zeros of the eigenvalue condition on the next Riemann sheet. We have 
made the assumption that the eigenvalues are real for arbitrary slab 
half -thicknesses. We have shown this to be true for thick slabs and it 
has been proved rigorously by others for the above-mentioned special 
cases. On the basis of our sample calculations, we conclude that if one 
is given the material properties Cj, and slab half -thickness, a, 

as well as the bare-slab eigenvalues corresponding to c 2 , <*2 811(1 8 
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then he can conclude whether or not the set ( s ») is empty and the 
maximum number of s n in { s n)‘ 

Perhaps the present results can serve as a guide for a rigorous 
Lehner-Wing type analysis of the reflected-slab problem. If the 
eigenvalues are all real, then one might be able to prove it in such an 
analysis of the present problem. 
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X. APPENDICES 


A. Summary of Elementary Solution Properties 

In this Appendix, the elementary solution properties derived by 
others (refs. 1, 4, 8, 9)t following the lead of Case (refs. 5, 7), are 
summarized. These solutions are obtained from Eqs. (2.1 6) - (2.18) of 
the text and are given by Eqs. (2.19) - (2.21). Such solutions are 
complete and orthogonal in the following sense. A function, say g(v<)> 
satisfying very weak restrictions (see for example ref. 7, Appendix G) 
for -l<a<(i<( 3 <l can be expanded as follows: 

1. Full Range; a = -1, 3=1 

8(t0 = [ a 4%j (M) + V-voj^V 2 ) 

+ f A j( v )<Pjsv^) dv ' (A.!) 


where the notation 6j(s) was defined by Eq. (2.25). The orthogonality 
relations used to determine the expansion coefficients in Eq. (A.l) are 

f ml Wjs V v >‘> > jsv(‘‘> dv * A jCv')vaj.(v)oj,(v) 

and for seS . . 

Ji 

/ x W Jsv .(n>9 tV0J (^ = °, 3 




’ 5 °j°j v 0j C js( ±v 0j)' 
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fi js( v Oj) = S n Js( z > | z =v 0J 


(A.3) 


2. Half Range; a = 0, p = 1 


8 W = a j%> )8 ;l (s) + X A J (v)<l ’js> )dv - (A-*) 


Here the orthogonality relations for seS^ are 


J 0 - A j(v , )w js (v')a5s(v')«j s (v'), 

/ 0 - vc j CT j v 0 jXj s ( - v 0J )<p_ voj ( V ) , 

/ 0 w Js(u)q> V0J (n)^ V0J (M)a^ = T fj^2i) a x js (±v 0J ), 

/ 0 W JS (H)<P Voj (^)^ s( .v)(,)d, = (^-) 2 vv 0J X Js (-v), 


f 0 V^sv'^jsf-v)^ 4 ** ■ Vtvoj + v , )Xj S (.v)<p JS (. v) (v’) 
and 

X ■ I c j <J j v > (A.: 


where 


W, (n) = — t w J J r 7 r, 0 < H < 1 , 

J 2fi Js (»)( Voj + H)X JS (-M) ~ ~ 


(A. 5a) 


(A. 6a) 
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?y>2 . n y ^ 

v* 0 


> o < m < i, 


with X 1q (z) given by 


V B) - ex p 


ni k 

In -*■ 

0 15 


4e( y ) dv l 
n;„(v) v-z| 


fij s (z ) 


(v 0j ' z2)n Js(“) x Js(- z ) 


c j a J 


— i * in»ua 

2fi js («) J 0 ( Vqj - n 2 )X Js f|i)(|i-z) 

The orthogonality relations for seSj e for the expansion (A. 4) are 

I 

S 0 A j( v )W* )dv = A j (v ' )W js (v ' )fi js (v,) S (v,) ' 

f 0 V^sv^^-v)^ * \ c J <T j v,X 0j S (- v ) < Pjs(-v)< v ' ) 


V^hsv^ 411 = I °j a j v ’ 


where 


°.iV 


W 1=! ()i.) = — r , 0 < (i < 1, 


Xj, (w) (w) 

E , 0 < H < 1, 

x ojs(t‘) «J s (m) 
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with Xqj s (z) given by 


X 0 js(*) 


exp 


L. r 1 ln 

4s(v) 

dvl 

| 2«i J 0 

Sk / A 

_0j s (v)_ 

V-Zj 


n js( z ) 


V” )X 0Js ( - Z > 


= l + 




/ 

J o 


(idjl 


2fi js (oo)^ 0 X 0js (-^)(M- Z ) 


(A. 8b) 


These half -range orthogonality relations and identities are obtained by 
extending the time -independent results of Kuscer, McCormick and 
Summerfield (ref. 15 ). 

A result, due to Kuscer and Zweifel (ref > 14), which we shall need 
to analytically continue solutions follows from Eqs. (A. 8). For a 
fixed value of z, Xj s (z) does not become Xqj s (z) as s crosses Cj. 
However, it follows from the middle expressions in (A. 8a) and (A. 8b) 
that 


( y 0j - z > x js( z > 


= 3 W !!) 


(A. 9) 


seS 


Jt 


seS 


je 


Following Kuscer and Zweifel (ref. 14) then, we define a function 
Xqj(z, s) which is continuous as s C ^ by Eq. (5.25). Such a 

function of the two complex variables z and s has the following 
analytical properties (ref. 14) . 

Fixed s: no singularity in z-plane cut along (0, 1); 

one simple zero at z = Vqj(s), Re(v Q j) > 0, 
only if seS. . 
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Fixed z: no singularity in the s-plane out along (-Oj, -ajd-Oj)); 

one simple zero at s » -Oj + c^o^z tanh“^(l/z) 
for Re(z) > 0. 


We note here that Xq^(z,s) is a nonvanishing analytic function of z 
and s for Re(z) < 0 and s^(-<Jj, -a^(l-c^)), the branch cut of 




82 


B. Derivation of 

In this Appendix, explicit forms of tjp+(x>H,s) are obtained. 
Following Bowden (ref. 1) we take, for medium j, the function 

Sjs( x >W x o) as 


- D j< x o)<P_v 0 J<• 1 > , 


V s > 


“(s+0 )(X-X Q )/V 
Cj(x 0 ,v)(p Jsv (n)e J 


dv, x < x. 


g Js (x,^x 0 ) = 


-(s+a J(x-x n )/v n , 

C J (x o)l>v 0 >) e J Je j( s ) 


T 1 / v / * -(s+Cj)(x-x 0 )/v 
+ J Cj(x 0 ,v)<p jsv (n)e dv, X > X Q . 

s. 

The expansion coefficients in Eq. (B.l) are to be determined so that 
gj S (x,n;xo) satisfies Eqs. (5-5) and (3 .6). That is, on putting the 
expansion (B.l) into Eq. (5 .6), we obtain in the limit x -» Xq 

f j(x 0 ,n)/4 - foj(x 0 )qp Voj (^) + D j(x 0 )«P. Voj (M)l 6j(s) 

+ f Cj(x 0 ,v)<Pj sv (n)dv. (B.! 

This is a full-range expansion (see Eq. (A.l)) of the function 
fj(xQ,|±)/n and use of the orthogonality relations (A. 2) gives the 


coefficients as 
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Upon vising Eq. (B.l), we obtain as 

i'2p(x,n,s) ■ jy 0 2 (x 0 )e^ +ff2 ' 02 dx 0 J » Vo2 (x,n,s)8 2 (s) 

^o] 

] W*> 


(s+<j 2 )x 0 /v 


-(s+a 2 )x 0 /v 02 


d x. 


H,s)dv 


|A,s)6 2 (s) 


+ f 0 [f* ° 2(x °' v)e 

- [/ x " Dg(x 0 ). 

pl r pa -(s+ff 2 )x 0 /v "| 

- J o LJ X C 2< x 0»-v)e ax 0 Jt 2( . v) (x,„s) d v. 


(B.6) 

The definite parity particular solution ^^(x, |i,s) is then obtained 
using Eq. (B.6) as 


V (x ’“' s) “{/I 


&2t( x 0> v 02^ e 


(s+a 2 )x 0 /v 02 


dx, 


r- x , (s+o 2 )x 0 /v 02j 

u J_ a ° 2 ±( x 0 ' v 02 )e * 0 , ,., 02 


)] *v 02 ( x '**' s ) 

®] W*’ 


M,s)>5 2 (s) 


>ir px 

J Cs±(x 0' 

u *a 


*)e 


(s+a 2 )x 0 /v 


dx, 


J *2v( x > 


H,s)dv 


± J'olf- * °2± (x O' v)e(S+02)XO/VaX o] +2(- V )( x '^ 


s)dv 


(B.T) 


That (B.7) is a solution of Eq. (2.13) for J » 2 can be seen by direct 
substitution as follows. The tj V (x,U,s) in Eq. (B.?) are solutions 
of Eq. (2.16), the homogeneous equation corresponding to Eq. (2.13). 
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However, their coefficients in Eq. (B.T) are functions of x so that 

some additional terms are obtained from the — operation. Thus we get 
r dx 

>* < c 2±(= t J ''02)%2 (ll) + C 2±( x »-''02) < I > .v 02 (t‘)]82(8) 




(B.8) 


which is an identity since according to Eq. (3*11) the are the 

full-range expansion coefficients of fg+fc,}!)/*!. 

To get ^p(x,|i,s) according to Eq. ( 5 . 7 ); we first note that 


dx 0 = 


• dx n + 


(medium) 

1 


0 * J * <*0 


X n -a 

• dx Q + / • dx 0 + 


»a ox 

» dXQ + / • dxQ + 


• dx Q , x < -a 


dxo. x > a. 




follow the same procedure as before and get ^ip+( x jK> s ) as 


(B-9) 
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t lp± (x,u,B) - 


and. 

\±U,H,s) = 


/ c l±( x 0 > v 01 ^ e 

¥ -cc 


(s+ffi)x 0 /v 01 


axJ» Voi (x,i*,B)6 1 (B) 


Jq[J X c &( x 0> v )^ to 0 Jt lv (x,n,s)dv 



-a 

Oy.(x 0 , -v 01 )e 


-(=+°x)x 0 /v 01 


dx, 


0 


a —a 

J C l±< x 0’ +V 01> e 

u -oo 


ax a - a 

/ -/ 0 1± (x 0 ,-v)( 

J o x 


+(s+a 1 )x 0 /v 01 


,] t - 


& 0 t. v (x,n,s)6 1 (s) 


01 


-(s+oJxq/V 


dx. 


n -a 

± / C li;( X 0 ' V ^ e 

^ mOO 


(S+(T 1 )X 0 /V 


*] 


dx rt |t-,/_ v j(x,4,s)dv, x. <-a 

(B.lOa) 


,-a 


•00 


C 1± (x 0 ,v 01 )e 


(s+Ol)x 0 /v 01 


dx. 


a -a 

+ / C 1± (x 0 ,-v 01 )e 

u „ x 

nl^A-a 

•Jo l/- Cl±(X °' V)e 


-(s+o^Xq/v -I 

^ 0 Uv^'X^,s) 8 1 (a) 


01 


(s+a 1 )x 0 /v 


dx. 


a -a 

T / C!±(x 0 , -v)e 

J -x 


-(s+d 1 )x 0 /v 


.] 


dx 0 t lv (x,n,s)dv 


" n-X 

± / C 1± (x 0 ,v 0 i)< 

Jr "«oo 


o]+-vqi 



(s+G^Xo/v 


dxoJ^_ Vm (x> s)&i(s ) 


x > a * 


(B.lOb) 
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Again, it is easily shown by direct substitution that Eqs. (B.lOa) and 
(B.lOb) are solutions of Eq. (2.15) for j - 1. We introduce the 
F functions of Eqs. ( 3 • 10 ) and by allowing x to take on negative and 
positive values, it follows that Eqs. (B.7)> (B.lOa) and (B.lOb) can 
be written as Eqs. (3*8) and (3*9) of the text. 

We also note here that the Cj+ coefficients of Eqs. (B.4) have 
the property 

C j±(“ x 0'" v ) = T C j±( x 0' v ) 

and 

c j±(- x 0>- v 0j) = T c j±( x 0> v 0j)> (B.ll) 

so that it then follows from Eqs. (3*10) that 

F2+(a,-0),s) = + F2f(a,0),s) 


and 


F+(-a,-n>,s) = + F±(-a,o>,s). 


(B.12) 
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C. Two -Media Full-Range Expansions and Orthogonality Relations 


In this Appendix, we first summarize some results of Erdmann 
(refs. 8, 9) and Kuscer, McCormick and Summerfield (ref. 15). Erdmann 
(refs. 8, 9) shows that a function, say h(ji), satisifying very weak 
restrictions for |i on the interval -1 £ [i < 1 can he expanded as 


h(n) = a 1 tPv 01 (^)5 1 (s) + 'b 2 <p_ Vo2 (n)5 2 (s) 


+ J A 1 (v)9 lsv (^)dv + J A 2 (v)<p 25V (jji)dv. (C.l) 


This is a two-media full-range expansion of the function h(p) and 
the expansion coefficients in it can he determined using orthogonality 
relations which are easily determined from the time -independent ones 
of Kuscer, McCormick, and Summerfield (ref. 15). For 
5-j_( s) = 5 2 (s) = 1 , that is s e A s 2i* these relations are 

f W s (|a)$ sv '(ji)d*x J' A(v)0> sv (n)dv = AtvOWstv'^v'^V'), 
f W s ([x)4> sv (p)9 Voi (|i)dii = 0, 

f ® sv ( i*) v_ Vo2 (**) a#* = 0, 

W s^Ko1^ )<P *v 0 2 (p ^ = °> 

f W s (p)% v (p)q>. Voi (4)<^ = vc(v)a(v)v 01 (v 02 - / v 01 )X s (-v 01 )<p. Vol (v), 
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f W s (|i)4 sv (n)q> Vo2 Ma|i = vc(v)o(v)v 02 (v 01 - v 02 )x s (v 02 )<p Vo 2 (v), 

4 

«!■ ? 
f l Ws(n) < T , v 01 <^> < Ptv 01 (^) 34 = - (v 01 ± VoaJXgiivo!), 

/ x W s(^) < f’.v 02 ^) < P±v 0 2^) a * i = ( C - g V ° g ) ( v 02 + v 0l' X s( ±v 02)> 

/ Ws(n)<P V01 (n)<Pv 02 (n)^ = - | VlWoi V 02 2 x s( v 02^ 


r« 

; -i 
and 

fw 

' -1 

where 


W s (^)q>. Vnp (M)9. V01 (^)dfi = | c l^ 1 c 2 ff 2 v 01 2 v Q2 X s( - v 01 ) > 


o(v),ff(v) - 


4 sv(n) ■ 


«4(v) = 


(C.2) 


X E (z) 


- X ls (z)X2 s (-z), 


t f ( v 02 + 

W 8 (v) = { 

Ufv™ - 


v 02 + v)Xg s (-v)W ls (v), v > 0 
(v 01 - v)X ls (v)W 2s (-v), v < 0. 


(0-5) 


All remaining quantities have been defined in Appendix A. 
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Rather than write out explicit orthogonality relations for the 
other three regions of the transform plane in the present notation, we 
introduce a function which is continuous as s From the results of 

Kuscer and Zweifel (ref. 1^) quoted in Appendix A, we see that one such 
function is given by Eq. (5*24) and can be written using Eq. (3*25) as 


Xq( z j s ) = < 


^ ( v 02 - z ) x 2s( z ) 

( v 01 * z )^ls( z ) > 

s € s n A Sgi 

(v 02 - z)]^ s (z) 

1 ’ 

s e S le A S 2 i 

A 

%>s( Z ) 

s e Sli A S 2e 

(v 0 i - z)x ls (z)' 


W z > 


8 


e ®ie ^ Spe* 


(C.k) 


In terms of this function, W s (v) can be written as 


c l*l v „ , , 

^)^- v ’ e) ’ v>0 


W s (v) = 


°2 a 2 v 


^2fi2 S (oo) Xq(v,s)' 


v < 0. 


(0.5) 


The function X s (z) is expressed as 


Xo(“Z,s) 


^ls( z ) 


( v 02 + z ) (v 0 l - z)J2 ls («>) 


■, Re(z) > 0 


* s ( z ) 


a 2s (z) 


( v 01 - z)Xo(z,s) (v 02 + z)fl2 B (»)' 


(C.6) 


Re(z) < 0. 
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In order to obtain a two -media expansion in the form (C.l), we 
generally have to switch some continuum solutions in one medium to those 
in the other. From the explicit form of 9j SV (l i ) (Eqs. (2.19 ) ) , we 


have that 


where 


0 l <, l < f2sv(^) - WlsvM = k s 8(u - v), 


•s s S (°1°X - =2»2) + - o 2 ) 


(0-7) 


We see that when the two media are the same, k s = 0. This quantity 
can be expressed in a number of different ways, and several that we 


shall use are 


{ c l a 1^2s( v ) “ c 2^1s( v ) 
c l cr l fi 2s( v 0l) 
-C2ff 2 ^i s (v 02 ). 


(C.8) 


The orthogonality relations (C.2) can now be written in terms of 
Xq(z,s) and k g as 

f W s (^)$ sv .(n)c4i J A(v)4>g V (|i)dv = A(v' , )W g (v , )^('V , )^s,( v< )^ 


W s (|i)^ s v(^)9v 01 (^)% * 0, 


w s M®sv(^)9.. v = 0, 

1 


J W S (w) < Pv 01 (^) < I'.v 02 (' i )^ - 
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vc(v)a(v)v 01 k s 


-1 01 4S22s( 00 )X 0 (-Voi ;S )(voi + v) 


f 1 w.(n) 




SVr“V'r S v\r/T V Q 2 


vc(v)g(v)v 02 k s Xo(-v 0 2;s) 
4 fi ls (a.)(v 0 2 - v) 


(|woi) 2 ^j2jloxIxo(-v 01 , s ) 


f w s(>‘) < i'v 01 (i*) < P ±Voi (M)ati = 


i w s^) ( P-v 02 ^) < P±v 02 (h)^ = 


G l cr l v 01 k s 

8fi 2s (o°)Xo(-v 0 i,s) ' 

" c 2 a 2 v 02 k s x o( ~ v 02> s ) 


% S M 


C 2°2 V 02 


\ 2 fl> 2s( v Qg) 1 

^2s(°°) Xq(-Vq 2 , s) 


f 1 W„(n)<f. ((.)? (M)dM .. c l g l v OlV02^Xo(-v og ,s) 

J -1 S ^ i+n ls («)(v nl - v 0 g) 


f w s(M) < f’. (n)ip. v (i»)au » r n — v g w 0lV ° gk 

02 01 _ ( 00 j ( Vr\r> — 


4S2 2s («,)(v 02 - v 01 ) X 0 (-v 01 ,s) 


These expressions appear more complicated than the corresponding ones 

in Eqs, (C.2); however, the orthogonality relations needed for all 

regions of the transform plane are given by Eqs, (C.9). That is, for 

s e Si e A S 2i , the proper orthogonality relations are the first, third, 

sixth, and eighth equations of (C.9) with Xo(z,s) given by Eqs. (C.4). 

We note here that Xq(z, s) always appears in Eqs. (C-9) with 

Re(z) <0. It follows then from Appendix A that for Re(z) < 0, 

Xq(z,s) is a nonvanishing analytic function of both z and s except 
for the branch cuts in the s -plane due to v Q1 (s) and v Q2 (s) . 


I 
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D* The Two -Media Full-Range Expansion for This Problem 

In this Appendix we show that application of the continuity 
condition, Eq. (2.15), results in a two-media full-range expansion of 
the type discussed in the last Appendix. For x = a, we readily obtain 
from Eqs. (5*1*0 and (2.15) upon using the explicit forms of tyj C+ . and 
tjp+ given by Eqs. (5*5), (5**0, (5*8) and (5*9) that 


0 = &2± [tv 02 (a,H,s) ± ^_v 02 (a,|i,s)]5 2 (s) 


+ S Q A 2±( v ) [| f 2v( a ^^ s ) ± t2(-v)( a ^^ s 3 dv 

? ^ a l±^v 01 (a^^s)6 1 (s) + A 1± (-v)\|r lv (a,|i,s)dvj 

+ F 2± (a,v 02 ,s)iJr Vo2 (a,p,s)5 2 (s) + J F 2fc (a,v,s)f 2v (a,p,s)dv 


- F2^(-a, v 01 , 



s)[t v oi(a,H,s) ± t.v 01 (a,kt,s)]6 1 (s) 
v,s) [t lv (a,p,s) ± t 1 (_ v )(a,|i,s)j dv. 


(D.l) 


We have indicated in Appendix C that according to Erdmann (ref. 8) the 
functions <Pv 01 (n)! <P-v 0£ >)> 0<V<1 and q> 2sv (n), 

-1 < v < 0 form a complete orthogonal set of basis functions for the 
expansion of h(|i), - 1 < p. < 1 for s e A S 2 i (see Eq. (C.l)). 

However, Eq. (D.l) also contains terms in which <p 2sv (|i), 0 S v S. 1, 

and cp lsv (u) , -1 < v < 0 appear. These continuum solutions must 

be replaced by corresponding continuum solutions for the other media* 
We use the relationship (C . T) to do this; that is, 
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and 


where 


■WtO - 

[sfsf VlsvM * ^ 6 < v - H(v) 


*U»(»0 

[c£3£ - 0^2 5 < v - “>] H (- v >' 

(D.2) 


f 1, v > 0 



H(v) = { 

(D-3) 


Lo, V < 0. 



When explicit forms of the elementary solutions and Eqs. (D.2) and 
(D. 5) are used in Eq. (D.l), we obtain the two-media full-range 
expansion 

* 

h(^) = fp^C-a/Vo^s) ± a3+Je"^ s+cr l) a / v 01 < P V()1 (jx)8 1 (s) 

¥ a 2± e ^ s+cr 2) a / v 02 <p_ v ^( n) 6g( s ) 

+ °C ± A i±(- v )j e ”^ s+ai ^ a ^ v 

- I^jpg+C^v^s) + A 2± (v)Je-( s+<J 2) a /v|9 lsv (^)dv 
± f° L [^ F;i±( ’ a, “ V,S ) e ' (S ' Wl)a/V ’ A 2t(" v ) e " (S+02)a/ ^ ]<P2sv(^)dv, 

(p t 4a) 


where h(p) is given by 
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h(|i) « ^(^(a^s) + A^Cfi)] e -( s+<J 2) 8 /^ 

c l a l 

± -*2_ H(-n)P 1± (-a,- l i,s)e-( 6+ °l) a /» 1 
C 2 <J 2 

+ [F2±(a,v 02 ,s) + a 2± "]e“( s ‘ Kr 2) a / v 02 cp Vo2 (|i)& 2 (s) 

T F ai .(-a,v 01 ,s)e (s+o l) a/v 01 . (D.lfb) 

The orthogonality relations, Eqs. (C.9)> can he used on the expansion 
(D.4) to obtain equations which determine the remaining unknown 
coefficients implicitly. However, it is convenient to introduce first 
the Ej + coefficients given by Eqs. (5*15) • We then have, after some 
algebra, the equations listed in the text as Eqs. (3-l6)-(3*23) • 



96 


E. Complex Representation of t s +(x,n) 

In this Appendix we outline how Eqs, (5,15) - (5*25) are extended 
to the complex plane (\x -* z) with s considered a parameter. As 
shown in references 2 and 3 , functions such as those introduced in 
Eq. (5*15) are extendable. The particular grouping of terms in Eqs, 
(5.16) - (5*25) indicates some integrals and residues which go together. 

The first functions to be considered are the Fj + (x,<o,s) functions 
given by Eqs. (3.IO). In Eqs. (3*20) - (3.23), these functions appear 
with Re(oo) >0 so we consider the functions Lj + (x,v,s) given by 
Eq. (5* 55)* When the explicit expressions of Eqs. (5*10) and (3*11 ) 
are used, we can show that for fj.f(xQ,fi) extendable n z without 
singularities in the finite z-plane then Lj+(x, v,s) can be extended 
to Lj+(x,z,s) given by Eqs, (3.32) and (3*55) • Now as z -> v e (0,1) 
it can be seen that the limiting values of Lj + , namely L* + (x,v,s) 
and lj + (x, v,s), are identical. Thus, L^ + does not inherit the branch 
cut of flj S (z) as one might be led to expect from Eq. (3*52) • There 
appear to be no other singularities of Lj + in the finite z-plane, 

Re(z) > 0 and Re(s) > Oy It follows from Eq. (5*52) that 


*i 

L j± (x ' v 0,j' s) = ~ c d^V OJi fi^ s (v ni )F^ ± (x,v ni ,s) 


2 


-(s+aj)x/v 0 j 


, S 6 


(E.l) 

In order to extend the functions Ij+(v) to the complex z-plane. 


we need the identity 
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Cg<*2 , ks f C l q l X 2s^ V ^ * C 2 g 2 X ls^ 


c 2 a 2 nj s (v)n £ s (v) C 1 0 1 c l a l c 2*2 fl£ s (v)q s (v) 


(E.2) 


which can be verified directly. We use this identity and find that 

Xj+(v), given by Eqs, (5*20) and (3.22), can be written respectively 

as Eqs. (5. 50) and (5.51) . The restriction Re(s)> -cr m on these 

equations comes from the fact that Lj± for both j = 1 and 2 occur 

in each I,. . More will be said about this restriction later. The 

j± 

contours C' are given in Figure 5* By letting z * vq 2 in Eq. (3,30) 
and z = Vq^ in Eq. (3*3l)> it can be seen that 


I 2±^ v 02 ,s ^ = J 2±^ v 02^ 


I l±( v 01' s ^ * J l+( v Ol)- 


(E.J) 


Thus, the inhomogeneous terms of Eqs, (3,16) - (3*19) are seen to be 
extendable and related as shown in Eq. (E.3), For z -> v 01 in Igf 
(see Eq, (3*30)) and z -4 Vq 2 * n l^i: (see Eq. (3.3l))> these 

functions might seem to be singular. However, it is seen upon examining 
the residues that this is not the case. Thus, the Ij ; j ; (z,s) appear 
to be analytic in the finite z -plane, Re(z) >0 and Re(s) > -0^. 

In Eq. (3»l6), we now let v -> z and for Re(s) > -a m and 
Re(z) > 0 in the finite z -plane, Eg±(z) is given by the inhomogeneous 
term 1 2 +(z), a term involving a 2 + if s s S 2 ^ and an integral over 
Eg+Xu), 0 < 4 < 1. A singularity occurs in the integrand when either 
fi 2s^) 0r fi 2s^ vanishes and this happens for s e c 2 » However, for 
this case, it is seen that we obtain from Eqs. (3.I6) and (3. IT) that 
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E 2 +( v 02) is related to a 2 +, It appears that Eg^(z) is analytic in 
the finite z-plane, Re(z) > 0, Re(s) > -a m and can be written as 
Eq. (3.28). We follow the same procedure with Eq. (3.17), and obtain 
an equation which is easily seen to be Eq. (3.28) evaluated at z = v Q2 ; 
that is, E 2 +( v 02> s ) an( i a 2± are related as 
1 . (s+(j 2 )a/v 02 

e 2± ( v 02> s ) ** 2 c 2 a 2 v 02%2s( v 02) a 2± e > s € ®2i* (EA) 

In a similar manner, we obtain from Eqs. (3*18) and (3*19) on letting 
v -> z and making use of Eq, (E.4), Eq. (3*29) of the text and again 
it follows that Ej4.(vqi>s) and a-j+ are related as 

1 , -(s+a 1 )a/v 01 

E l±( v 01> s ) = 2 c l a l v 01 fi ls ( V 01 ) a l± e > s e S li* ( E *5) 

It also appears that E 1 +(z,s) is analytic in the finite z -plane, 

Re(z) > 0 and Re(s ) > — cr m . 

The solutions tj C ± and 1^4. can now be written in terms of the 
E ,.j. as shown in Eqs. (3*3 6) - (3*39) Of the text. 

tJ““ 
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F. Investigation of the Associated Eigenvalue Problem 

We consider in this Appendix, the associated eigenvalue problem; 
that is, the problem for which f(x,u) so. The inhomogeneous terms 

(given by Eqs. (5*50) - (3*32)) can be seen to be identically zero 
when f (x, }i) is zero everywhere. Solutions for I 1± =0 will be 
denoted with a bar, i.e., Ej + . The unknown expansion coefficients for 
the eigenvalue problem, Ej ± , are given by Eqs. (5.28) and ( 5 . 29 ) vith 
Ij-f =0. It is seen from such equations that Ej+ can be determined 
only to within an arbitrary factor independent of z and that 
depends on Eg*.. Furthermore, the original normal-mode expansion 
coefficients for the eigenvalue problem are given by E^ + ((i,s), 

0 < H < 1 , j = 1 , 2 ; E^v^s), s e S^ and ^(v^s), s e S 2i . 

Therefore we must examine solutions of such equations as a function of 
the transform variable s for z -> p with the contour C* collapsed 
onto the branch cut ( 0 ,l) due to fl 2s (z' ), and for z = when 

s e S . . . This will be done for all s in some right -half s -plane and 

J 1 

it is convenient to divide the plane into three regions: s e Sg e , 

s e S 2i and s e C 2 « 

When s e S 2e , fi 2s (z*) does not vanish within C* so that Eq. 
( 5 . 28 ) with Ig*. = 0 can be written as 


E 2 f(|i,s) = ± K g (|i, v)Egf (v,s)dv, s e Sg e , 0 < U < 1, 


(F.l) 


where 


K s (h,v) 


k r, fi 2s(°°) Xq( -u, s )X q( -v, s ) v -2(s+a 2 )a/v 

- — rr e , 0 < |i> v < 1, 

2 fi^ s (v)n- s (v)(v + 4) ■ " 

t* ” v 02* 


(F.2) 


I 
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When s e S 2i , n 2s (z') vanishes inside c' hut not on the real 
interval (0,1). As C' is collapsed onto (0,1), a residue term appears 
so that Eq. (5*28) with - 0 takes the form 


_ , , k s fi 2s <“> x o (* | a ' s ) X 0 ( - v 02 ’ s) - 2 ( E+ ° 2 ) a / v 02 - 

Eo+(u,s) = ± ■ ■■ T T " ' ' e E' 

c 2 cr 2 fi .lsM fi 2s(v 02 )(vo2+fi) 

± jT K !? ((i,v)E 2± (v,s)dv, s € S 2i , 0 < W < 1, 


®2±( V 02 ,S ^ 


(F.3) 


for z where K g (n, v) is given by Eq. (F.2). However, Eq. (3*28) 
with I2+ = 0 must also hold at z = Vq 2 and this gives an additional 
constraint on solutions of (F.3)> namely 


. x k s » 2s <”> x o(- v 02’ s) - 2 ( S+Q 2 )a / V 02 - , 

E 2±<V’ s > ’ ± a** 2v 0g ^ s (v 0g ) 6 E2+ - (V °2> s) 


M v 02' v )®2±( v ' s ) dv J s e s 2i* 


(FA) 


When s e > 2 , the curve separating S 2 ^ and S 2e , ( v ) ft 2s ( v ) 
vanishes for some v on the interval (0,1). That is, Vq 2 is real 
and lies on (0,1). Setting Vq 2 = q, we can put Eq. (3*28) with 


=0 in the form 


, k c VoA 00 ) X 0 (-n,s)X 0 (-T),s) -2(s+a 2 )a/r) _ 

E (|i,s) = ± ± — e E 2 ±(t),s) 

2 ~ 2 c 2 a 2 flis(°°) ^ 2 s (t))(ti+h) 


ks fl 2s(°°^ 
2 




1 E 2 ± (v,s)Xq(-v,s) -2(s+cr 2 )a/v 


2s( v ) n 2s ( v > 


s e C 2 , 0 < fi, n < l. (F*5) 
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Note that this equation would be obtained from either Eq. (F.l) or Eq. 
(F . 3 ) for s -» Cg from s e S 2 e or s e Sgi respectively. 

For arbitrary complex values of s, the kernel K s (n, v) which 
appears in Eqs. (F.l) and (F. 3 ) is not symmetric since 


K s (n,v) + [*.(**>] , 3to(s) i 0, (F. 6 ) 

where * denotes complex conjugation. Note, however, that when 
Ira(s) = 0, the unknown functions Eg^hjs) can be redefined so that 
a symmetric kernel is obtained. Solutions of Eqs. (F.l) and (F. 3 ) 
depend on the behavior of K g (n, v) and we shall look at a quantity 
B 2 (s) given by 

b2(s> . r 1 / 1 

^0 J o 


2 

K s (n,v) d|idv. 


(P.7) 


To do this, we introduce the nondimens ional parameters £, a^, <?p 

and A given by Eq. (6.1) with £ = a + 1(3. Note that a, p, 0 p, Up 
and A are real while tfp and A are nonnegative. In terms of these 


quantities, we have 



(F.8) 









A 
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To make estimates of the function |Xq(-|i, s >r , we use an integral 
representation of the single-medium X-function given by Kuscer and 
Zweifel (ref. 14-), namely, 


f l r ico ) A «i 

x 0j(-^ s ) = ex P IpZT / ln ■ ^ | S / v — — 


(F.9) 


Upon letting z' = iy and using flj s (z’ ) = ftj S (-z')> we see that Eq. 
(F.9) becomes 


x 0j ( ' M,s) = exp ^/ c 


n “ ta n .i^ (ly) l 

o fi js(°°) y^ + n j 9 


(F.10) 


which is seen to be real for s real. In terms of the quantities of 
Eq. (6.1) we find 

I 2 2 -v 

i v , J 2 Jan r , \ [ a *°D-' J Re(yO + e> dy \ 

x 01 (-n, s )| = exp <— / In \ — 7 , g ng y g > 

1 J 0 V (a+c D -c R ) + p^ y 2 + H^J 


x Op(-^s) = exp 


2 m- r°° C a ~s(y )] + P 2 dj 


3t V (a-1) 2 + P 2 y 2 + M 2 j 


(F.ll) 


where 


g(y) - y tan l/y. 


(F.12) 


It follows from Eq. (F.12) that 0 < g(y) <1 for 0 ^ y < °° and that 
it is a monotonic increasing function of y on this interval. Further- 


more, since 


2E / *y 
n J o y 2 + n 2 


= i, 


(F.13) 
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i l 2 

the following bounds (perhaps rather loose) are obtained for |^o j ( "l 1 * s ) | : 


2 ■ e 2 




and 


2 I (o+o f + 0 2 

1 < X 01 (.|i,s) < \|- s £ 5 

1 01 1 N (-WOb-^r) 2 + 


(F.14) 


for a > max (l, -Ojj+a R ), and 


\) (a-l)i ; p - 


2 , q2 

a + p 


(a-l) 2 + P 2 


and 


Y 


P 2 ' I * x| 2 ^ /- \ I ( a+a D) 2 + 0 2 

< |X 0 i(-4,s)| < max I 1, 


2 . «2 


(o+a D -a R ) + p 


2 . «2 


(a+Op-aRr + 3 


(F.15) 


for a < max (l, -CTp+cr R ) and p ^ 0. We note that in the £ -plane, the 

points £ = 1 and £ = -<jj) + 0 R are the right ends of the branch cuts 

of Vq 2 and v Q1 respectively and these cuts lie on the real 

£-axis. The left ends are at £ = 0 and £ = -c D , respectively. 

| fl | s ( v ) I 2 

The functions — r=~ — I are easily found to be 

I c 2<*2 I 


«| s ( v ) 


c 2°2 


= (<x-v tanh"^v) 2 + ^p ± , 


(F.16) 


* 

Recall from Appendix A that the curve C 2 (see Figure 2) is given by 


a' = tanh“l . 
x x 


(F.17) 
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The parametric form of this equation is 


P’l = a' = v tanh - ^ v, 0 < v < 1. 


(P.18) 


We see from Eq. (P.l6) that 


c 2 a 2 


are the squares of the distance 


(in the £ -plane) from the point (<x,p) to the points (a'(v), + p’(v)) 
respectively which lie on the curve C 2 . Since these functions appear* 
in the denominator of |K s (|i, v)| the integral (F.7) will not be 
bounded when a and p * are related as in Eq. (P.17). We define 
I\ni n ( a ,p) as the minimum distance from the point (a,p) to the curve 


C 2 for 0 < v < 1. That is. 




c 2 or 2 


n | s ( v) 

c 2 o 2 


, 0 < v < 1 


(F.19) 


and rfoin(<x,p) / 0 for (a,p) £ C 2< 
(F.16) and (F.19) that 


Therefore , we have from Eqs . 


-1 <-r i , 0 < V < 1 . 

«2s(v)| 2 | S (v) r D min^P) 

C2 ff 2 I I 0202 I 


(F.20) 


Analytical bounds for the above function are not as easy to get. For 


0 = 0, fit (v) = (fi| s (vi] and 


«5 a (v) 


a 2 g (“’ V )’ 


(F.21) 


where g vj has been investigated and tabulated by Case, de Hoffmann, 

and Placzek (ref. 6). They show that g (— , v) I occurs at v = 0 

' V a /I max 


occurs at v = 0 


3.05 


for a < jt 2 /8 whereas for a > jt 2 /8 it occurs for v between 0 and 

1 , For a very large they have g max -> 4a 2 /jt 2 . The present geometric 

interpretation is consistent with all of these characteristics. The 

radius of curvature of the curve C 2 given by Eq. (F. 17 ) is ji 2 /8 at 

(a', (3') = (0,0). For a‘ very large, (3'-» jt/2 so that the minimum 

squared distance from (a, 0 ) to (a',p' ) approaches it 2 A, in agreement 

with Eq. (F.21) and g -» 4a 2 /jt 2 . 

max 

-4oA/v 

Note that the exponential factor e in Eq. (F. 8 ) requires 

a > 0 in order for B 2 (a,|3) to be bounded since both v and A are 
nonnegative. On using estimates (F.l4), (F. 15 ) and (F.20), we obtain 
from Eq. (F. 7 ) a bound for B 2 (a,|3) which we denote as B 2 __(a,J 3 ): 


-4aA 

e 


([a(a B -l).a D ] g + p 2 (g R -l) £ } (a 2 +P 2 ) 


(a+a D -£y R ) 2 + 


a > max(l, -Ojj+a^) 


e -4oA / a 2 + p 2 > 

— r max — — — — — 1 

V(a-l) 2 + P 2 / 


(a-l)* + p‘ 


^(gp-lj-gj 2 + p 2 (g R -i)^ Qa-1) 2 + B 2 J 


■, P / 0, 


0 < a < max(l,- 0 p+ 0 R ) . 


(F.22) 


io 6 


It can be seen that depends not only on a and (3 but also 
the nondlraensional material parameters and a R as well as the 
slab thickness parameter A. The estimate (F.22) for is not 
bounded for the following regions in the s -plane; 


and 


Re(s) < -Og (a < 0), 

S 6 C 2 

s e branch cuts of Vq^(s) \J v@ 2 (s). 


(P-25) 


These regions must be handled separately. Even for the general case, 
where s does not belong to any of the regions (F.25), it appears 
difficult to say whether or not the eigenvalue problem has nontrivial 
solutions. We suspect that it has only trivial solutions for such 
regions since that is the result which has been found for certain 
special cases by others. Lehner and Wing (refs. 1 6, 17) have shown this 
for the bare slab, while Lehner (ref. 15) and Hintz (ref. 10) have 
obtained this result for the slab surrounded by pure absorbers. We 
can show that this result is also obtained for the special case A -> 
that is, a thick slab. 

Since the slab thickness parameter A appears only in the 
exponential term of Eq. (F.22), it is seen that Bjnax^P) can be made 
as small as one likes as A -» « if s does not belong to any of the 
regions given in (F.23). For |®max( a '>P)| < 1, the Neumann series 

solution of the inhomogeneous integral Eq. (F.3) converges to a unique 

* 

solution. (See ref. 19, for example.) Fredholm's Alternative Theorem 
(ref. 19) then guarantees that the corresponding homogeneous equation, 
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namely (F.l), has only the trivial solution, Thus, for s e S 2e the 
eigenvalue problem has only the trivial solution as A ». When 

s e S 2; j_, the unique Neumann series solution of Eq. (F. 5) must satisfy 
the additional constraint, Eq. (F.4) . Using the condition fi 2s ( v o 2 ) a 0, 


we obtain 


-2(s+02)a/v o2 i0( s 

e = p(s)e 


r. 


(F.24) 


where 


P 2 (s) 


pe(v 02 )-l] + Im 2 (v 02 ) 
[Re(v )+l] 2 + Im 2 (v 02 ) 


and 


0(s) = tan"^- 


Im(v 02 ) 


^ e ( v 02^ ” ^ 


- tan 


-1 


Im(v 02 ) 


Re(vo 2 ) + 1 


(P.25) 


Now since Re(v 02 ) > 0, we have 


r p<s) ] °» Re(v o2 ) ^ °- 


(F. 26 ) 


^ ' J A-»00 

Therefore, the Neumann series solution is seen to converge to zero as 
A -> «> when Re(vQ 2 ) / 0. Note that Re(vQ2) = 0 is the branch cut 
of v 02( s ^ which is one of the regions given by (F.23) which we must 
consider separately. When s e C 2 , v - tj, 0 < tj < 1 so that p 2 ( s) 
of Eq. (F.25) becomes 

p2(s) = (nTl) - 1 (F-27) 

and p = 1 occurs only at s = -o 2 (that is, (a,0) = (0,0)). 


~V 
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We use Eq,. (P.27) in Eq. (F.5) and on talcing the limit A » we find 

¥. 


that Eg4.(y,s) 0 for s e Cg, s ^ -Cg. 

Summarizing the results then for A -» «, we find that the eigenvalue 
problem has only the trivial solution for Re(s)> -<jg unless e 
belongs to either the branch cut of v 01 (s) or v Q2 ( s), In order to 
see what happens on these cuts, we must write Eqs. (3.28) and (3.29) 
with Ij ± =0 in terms of the Xj s (-z) functions rather than the 
Xq(-z,s) function. This will be done in the next Appendix. When A 
is not large, one has for the problems of Lehner and Wing (refs. 16, 

17), Lehner (ref. 15 ) and Hintz (ref. 10) that if the eigenvalue problem 
has nontrivial discrete solutions, they occur on the real s-axls. For 
the bare slab, Bowden (refs. 1, 4) has shown that these solutions lie 
on the branch cut of Vq(s). In view of these results, it is assumed 
that the eigenvalue problem has nontrivial solutions for Re(s) > -Cg 
only if s belongs to either the branch cut of Vq^(s) or Vog(s). 
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0. Solution of the Associated Eigenvalue Problem for s c Sgj 

In this appendix we look at solutions of Eqs. (3,28) and (3*29) 
with Ij ± s 0 for s on the branch cuts of v Q j(s) . It is convenient 
to use coefficients related to the original expansion coefficients 
Aj+(v) and (the bar Indicates that we are considering the 

associated eigenvalue problem). Recall that the Ej ± are related to 
such coefficients by Eqs. (3*15); (E.4) and (E.5)- We also noted in 
Appendix P that the coefficients can be determined only to within an 
arbitrary factor Independent of v. Following Bowden (refs. 1, 4), we 
introduce coefficients Bj+ as 

A j±( v ) * a st B j±( v )> a l± * a 2fc t 3± * (®*1) 

p 

The estimate B mf>v (a^p)^ Eq. (P.22) of Appendix F, was not bounded on 
the branch cuts of Vqj(s). In that estimate we used the Xq(-z,s) 
function so that the behavior for s inside, on, and outside the curve 
Cq could be seen. To investigate what happens on the branch cuts of 
Vqj(s), we should use the Xj s (z) functions (Appendix A) which do not 
inherit the branch cuts of Vqj . Also, when Vq^(s) becomes pure 
imaginary (that is, on its branch cut), we cannot Include its contri- 
bution (the pole at z ' * v Q1 ) in the integral over the contour C 1 of 
equation (3*29)* We note again that the material properties Cj and 
<jj determine where on the real s-axis the branch cuts of Vqj(s) lie. 
The only restriction which has been made is c^ < 1 and this alone 

does not specify how the cuts overlap. It does, however, guarantee 

* 

that the branch cut of v Q1 (s) lies entirely to the left of s * 0. 
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We consider s s Sg^ A first. When the relationships (G.l), 
(5*15); (E.4) and (E.5) are used in Eqs. (F.3) and (F.4), we obtain 
after some algebra and use of the X-identities of Appendix A, 
equations for Bg+(n) and the additional constraint, namely Eqs. (4.2) 
and (4.6) of the text. Recall that Eqs. (F.3) and (P.4) were obtained 
from Eq. (3*28) with Ig± s 0. Equations for B-j^-n) and b^ are 
obtained in a similar manner from Eq. (3*2 9) with I 3± s0 when the 
contour C* is collapsed onto the interval (0,1) of the branch cut of 
^2s( z ')‘ These equations are given as Eqs. (4.3) and (4.4) of the text. 
The normal-mode expansion of the solution of the associated eigenvalue 
problem is given in terms of the Bj+ coefficients by Eq. (4.1). We 
note that Eq. (4.6) is the exact eigenvalue condition since all material 
properties have been assumed known. It determines the values of 
s, { s n}* ^ or which the eigenvalue problem has nontrivial solutions. 

When s belongs to tile branch cut of Vq^, Eq. (4.6) takes on different 
values above and below the Vq^ cut. Therefore, it is concluded that 
the eigenvalue problem has only the trivial solution on the branch cut 
of Vq 3 _. On that portion of the branch cut of Vq 2 which is not a3 ( - 1 
part of the Vq^ cut, Eqs. (4. 2) -(4.6) require that the limiting values 
of the coefficients above and below the vq 2 cut be related as 

[Bj ± (n)] + - ±[b j± (m3 ‘ 

and 

[b}+] * -[pi±] > Re(vQg) * Dn(v 01 ) m 0, (G*2) 










Ill 


that is, where s is real and is given by 

max £-(jg, -cri(l - < s < -a 2 (l - 09) • It then follows from Eqs. (4.1) 

and (G.2) that the limiting values of \jr s+ (x,|i) for the same region are 
given by Eq. (4.21). From Bowden's results (refs. 1, 4) for the bare 
slab, it is expected that the eigenvalue problem has nontrivial solutions 
only at isolated points, which lie on the branch cut of Vq 2 but 

not on the branch cut of v Q ^. 

In the limit c 2 cr 2 a 00 which was discussed in Appendix F, we see 
that Eq. (4.2) gives B 2+ (n) 0 while Eq. (4.6), the eigenvalue 

condition, becomes 

„ X 2s (-v 02 ) e - (s+o 2) a / v 02 X 2s (v 02 ) e ( s +° 2 )a/v 02 

U = ' " 11 ■ ’ ’■ - ■ + 7 ■ i — , CpOpa — > 00, 

x ls(“ v 02) V 01 + v 02 x ls( v G2) V 01 - v 02 


Re(v 02 ) * lm(v 01 ) = 0. 


(G.3) 

Equation (G. 3) is the "thick-slab" eigenvalue condition and for the 
region of the s -plane where it is valid, it can be seen that we have an 
even eigenvalue s n if 

’X 2s <v 02 ) e ( s “ +o 2) a / v 02‘ 


Im 


fX Se (v 

l x ls( v - 


02) ( v 01 


iS)a/v 0 2> 

-J = ° 


and an odd eigenvalue s n if 


Re 


r x 2s (v 0g ) .KggWgg i 

\X l3 (v 02 ) (v 01 - v 02 ) J 


= 0, c 2 cJ 2 a ~ >00 ' 


(0.4) 


We note that Eq. (G.3) has the same form as the zero-order approximation 
of the critical condition given by Case and Zweifel (ref. 7) , except 


m 
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here we have both even and odd parity solutions. Numerical solutions of 
the eigenvalue conditions will be discussed in Appendix J. 

In the region of the s -plane s € S^ e A we are specifically 
interested in the solution on the branch cut of VQgCs) which lies to 
the left of s = -crj. That is, for s real and -n 2 < s — *' or 

such values of s the solution (4.1) outside the slab is not bounded 

as x since 

t lv (x,|i,s) =<P l 6 V (n) e -( s+ 0 l) x / v , 0<v<l. (0.5) 

In addition, the restriction Re(s) > -a m on both inhomogeneous terms 
I (see Eqs. (5*50) and ( 3 . 31)) also indicates that we cannot deform 
the inversion contour to the left of Re(s) = -a^ in general. However, 
when one is looking for the solution inside the slab, |x| < a, perhaps 
the inversion contour can be deformed to the left of Re(‘') » -cr^ for 
special values of material properties and/or initial data. For 
s e S 2 i f\ S^e, expansion coefficients for the solution inside the 
slab are obtained as Eqs. (4.7) and (4.8). We note that Eq. (4.8) is 

exactly Eq. (4.6) with Xq 1 s (z) replacing (vq^ - z)X ls (z). Recall 

from Eq. (A. 9) that these are the X-functioris which are continuous as 
s —> C]_. Under the same replacement of Xq 1 s (z) with (v 01 - z)X ls (z) , 
Eq. (4.7) reduces to the equation from which Eq. (4.2) was obtained. 
Equation (4.8), which corresponds to the eigenvalue condition Eq. (4.6), 
determines the pseudo-eigenvalues. T t is, the values of s, 

-<j 2 < s < -a where s ) has nontrivial solutions. 
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H. Form of \jr S 4 .(x,ti) on the Branch Cuts of Vqj(s) 

In this appendix, we show how the transformed solution \|r s ±(x,|i) 
is put in a form where one can see how it "behaves on the branch cuts of 
Vqj(s). We expect that \|r g+ inherits the branch cut of v Q1 (s) since 
only one of the two discrete modes appears for jx j > a. Such branch 

VV 

cuts appeared in the half-space albedo problem solved by Kuscer and 
ZWeifel (ref. l4) as well as the two dissimilar semi-infinite media 
problem solved by Erdmann (refs. 8, 9)* We also expect that the branch 
cut of Vq 2 (s) is not inherited by \|r g + but Instead one should find 
poles at s = s n , the place where the associated eigenvalue problem has 
nontrivial solutions. This is what Lehner and Wing (refs. 1 6 , 17) and 
Bowden (refs. 1, 4) found for the bare slab. 

It is not obvious upon looking at the equations of section III which 
determine the expansion coefficients implicitly how we should group 
terms to show what we expect. We start by looking at \{< >> -j - (x,{i,s) . 

From Eqs. (5*3) and (5*8) of the text we have that 


t2±( x '^ s ) = a ; 


at [tv 02 ( x ^ s ) ± '|r- V02 (*,H,s)l 

+ A 2±( y )j^2v( x ^^ s ) * ^2(-v)( x >^ s )J dv 

+ F £± (x,v 02 ,s)\lr Vo2 (x,n,s) 

± ( -x, v 02 , s)t_v 02 (x, |i, s) 

+ f Q |Fah(x,v,s)\Jf 2v (x,ji,s) ± F 2± (-x,v,s)t 2 ^ v ^(x,p,s) 


dv. 


(H.l) 


40 
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Note that this equation can, he obtained of course from Eqs. (3*36) and 
(3>36)" It is readily shown from the definition of the functions, 

Eqs. (3*10), and the properties of the Cj + , Eqs. (B.ll), that 

F2+( x > v 02> s ) = Fp+( a ; v npj s ) - F2+(-x, -V q 2, s) 

and 

Fg±(x, -v Q2 , s) ^ ^2±( a f ” v 02> s ) ± F 2 ±("X,v 02 ,s). (H.2) 


It follows from Eqs. (H.2) then that two coefficients in Eq. (H.l) can 
be written as 

F 2 f(x,v 02 ,s) * | |f s± (x,v 02 ,s) ± F^-x^Vq^s^ + | F 2± (a,v 02 ,s) 


and 


±p a;(- x < v 0 S' s ) =|h(- - v 02’ s ) ± F S±(‘ x ' v 02' 6 )j ± |F2±( a > v 02< s )> 

(H.3) 

where we have used Eq. (B. 12) to replace F^. (a, -v 02 , s) . Equation (H.l) 
becomes then 

t2±(x,h,s) = j^a^j. + | s)J^ Vq 2 (x, h, s) ± l_ V()2 (x,n,s)J 

+ J o A 2±( v )|t2v( x ^> s ) ± t 2 („ v )(x,^s)jdv 
+ |^ F 2±( X ' V 02' S ) * F 2 ±(- x ,-v 02 ,s)JiJrv 02 ( x r^s) 

+ |^F 2± (x,-v 02 ,s) ± F 2± (-x,v 02 ,s)j\|r. Vo2 (x,|i,s) 

+ f Q p2fc( x ^ v / s )^2v( x ^^ s ) * F 2 ± (-x,v,s)\|t 2 ^ v j(x,^,s)| dv. 

» ' c ■ 

(H.4) 
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When s e branch cut of v Q2 (s), we have v 02 = 1 | v 02 1 for ^to(s) = 0“ 
and v Q2 = -i | v Q2 1 for Im(s) * 0 + . Therefore, on going from below 
the Vq 2 cut to above it, we see that the third and fourth terms in the 
RHS of Eq. (H.4) simply interchange while those containing P 2 +(x,v,s) 
and Fp 4 .(-x,v,s) are unaffected since these functions do not depend on 
v 02 . The coefficient of [a^ + | F*(a,v 02 ,s) 
for odd-parity solutions and we do not yet know how A 2+ (v) behaves, 

By comparing Eq. (H.4) with Eq. (4.1) for J x | < a, we suspect that 
£a 2+ + q F 2±( a > v 02> s )j ^ tie coefficient which excites the associated 
eigen -solution \Jr s 4 .(x,n) . This is the information we needed to see how 
to group terms in the implicit equations for the expansion coefficients. 

Now we look at the equations which determine the expansion 
coefficients. We obtain from Eq. (5*17) upon using the X-identities, 
the definition of the hj functions (Eqs. (4.5)), and the relationship 
between the E 1+ and the original expansion coefficients that 


changes sign however 


s |^ a 2± + |F£fc(a,v 02 ,s)j 
+ |P 2± ( a ,v 02 ,s) £ 


02) ± h 2(~ 
v 02 V 01 

h 2 ( v Q2) h2( _v 02) 


V 01 * v 02 V 01 - v 02 


- v 02) 1 
" v 02 J 

] 


[A*(,) + ( a, u, s )j h 2 ( n) 


01 


± f F l±(“ a >^ s ) h l(l 1 ) 


M + Vqi 

0 ~ |i 2 - v 02 2 


c4i 


± F l±( ~ a ’ v 01> s ) h l^ v 01^ 


2v 


01 


2 2 
V 01 * v 02 


(H.5) 
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Following the same procedure with Eq. (3.16), we obtain after several 
partial fractionings and use of Eq. (H. 5 ) 

= ||F 1± (-a, v , S ) e ( ^ )a A 

liii rt t u\ * 1 _ ^ % 1^" Irt _ l « ■ \ 


k s flgsH VQ2 - v M v ) 
^ %bW v oi^ “ 


h 2( v 02) , h 2^" V 02^" 


+ | F 2± (a,v 02 ,s) 


if 1 _ t .1 f h 2( v 02) , h 2^’ v ( 

x < K + 5 F ^ (a ’ v ° 2 ’ S >J LtT^ 1 — 

+ i P 2+ (a,v 02 ,s) Nl°il T h 2 ( - V 02 ) l 
g 2+' 02’ |_ v + v 02 + V - V Q2 j 

*f Q f A 2hM + F 2t < a ,M^,s)J hg(n) 

r 1 „ , , % v 01 2 - ^ 2*28! 

dn ■ v 2 _ „2 Cr>a, 


v 01 - H 2<P 2 sv (*i) 

v 02 2 - n 2 °2®2 V 


dv> . 


From Eq. (3*19)> we get 


?h 1 (-v 01 )a 1± =^a Sfc + |F sfe (a,v 02> s)|^ 2 (v 02 ) ± h 2 (-v 02 )J 

+ |%(a,v 02 , 8 )[ h2 (v 02 ) Th2 (-v 02 )] 

+ (^oa 2 - v°! 2) / o [ A 2 ±(^) + F 2± (a,fi,s)j h 2 (p)- 

(v 02 2 - v Q1 2 ) f F^-a^s^Cn) - 

'-'o H - v 02 ^ 


Fi±( -a, v Q1 , s) -v 01 ) ± h^v 01 )J 


(h-t) 
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Finally, we have from Eq. (5*18) that 


*„<-> » s -.i* 1 '- 1 * - gs is} 

*{[*»* l'a<*'*oa.‘l][SSi* Srf}] 


.1™ /_ .. ^ rM v O£> _ h 2 (-V02)l 
+ 2 P 2±(^v 02 ,s) [7—^ + v + v 0£ J 

*f 0 [ A at (,l) + F a±< a ^^ s >j h s<^) ^ 

t / 1 P 3± (-a,^B)h 1 (^) Sj-ij . g nrvl 


Vpi 2 - dfi 

v 2 2 |i + v 

v 02 “ ^ 


+ L(-a,v,s) - ^r a (a,v,Oe (l 'l- I 2 l8/ j . (l 

^ X J* 

We now follow the procedure of Bowden and Williams (ref. 4) and 
write the expansion coefficients Aj+(|i) and a-j* in the form 

Aj+(n) = |ap4- + ijjr ^2t( a>v 02 , ®)J^ r j+Cm-) + ®j±^) 


0*0 1 + 


|.F 2± (a,v 02 ,s)ja , a± + h 3± (^). 


When Eqs. (H.9) are used in Eqs. (H.6)-(H.8) it follows that 


A' j±(m-) = a 'i+ = ^2+, (H.10) 

where Bj + (|i) and 132+ are the expansion coefficients of the 

associated eigenvalue problem given by Eqs. (4. 2) -(4. If). The 

< 

coefficients B«+(y) and b-j+ are found to be given by Eqs. (4. 10) -(4. 12) . 
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The coefficient is obtained from Eg,. (H.5) as 

fa^ + | P 2± (a,v 02 ,s)l V °^ + r B 2± (n)h 2 (n) — ^ ] 

L 2 u JL V 01 + v 02 V 01 “ v 02 Jo 2 * + V 01J 


h g( v 02) h g(- v 02) 

V 01 + v 02 * V 01 " v 02 


* ■ j F a± {il ' U 02' L ' [: 

[■ 2+(M-) + F£f(a,(.i,s)J h 2 (n) - ^ 

LX 4* Va*i 

+ / P3i(-a,n, s)h x (yi) <4* 

J 0 M - V n o 2 


] 


01 


T Fi+(-a, v oi^ s )hi( v oi) 


2v, 


01 


v 01 2 - v 02 s 


(h. 11) 


It can be seen from Eg. (4.6) that the coefficient of 
£a- 2 ± + ^ F£ 4 .(a,VQ 2 ,s)J i n Eg. (H.ll) is the eigenvalue condition, and 
it will be zero at the places where the associated eigenvalue problem 
has nontrivial solutions. Equation (H.ll) appears in the text as 
Eg. (4.15). The solutions \|r,..(x,|i,s) can now be written as 

j-i- 

Egs. (4.16) and (4.17) of the text. 
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I. Behavior of ty B+ (x,n) on Inversion Contours 

In this appendix, several points concerning the behavior of 
\|r s+ (x,n) on the integration contour of the inverse Laplace transforma- 
tion and some portions of related deformed contours are discussed. 

First, we look at the behavior of \Jr 6+ (x,n) as | s | -» co with 
Re(s) = 7, a large finite positive number. It will be seen that 
\jr s+ (x,|i) is not necessarily 0^. Such parts of \|r s+ (x,M.) are 
inverted separately and the resulting solutions are shown to satisfy the 
uncollided transport equation. Then we consider how \Jr s+ (x,|i) minus 
the uncollided term, ''l' u +(x,fi>s), can be deformed around the poles and 
branch cut of \Jr s± which were discussed in section IV. 

We are interested in the behavior of \|r s+ on the contour 
Re(s) * y as | s | where y is finite. For such cases, 

s e S^ e A S 2 e and the solutions can be seen from Eqs. (3*3)> 

(3. if), (3-8) and (3*9) to be 

tsh(x,h,s) £ Agj^v) + F 2± (x,v,s)Je"^ s ‘ Hy 2) x / v q) gsv (| I ) dv 

± f Q f A 2fc( v ) + F^(-x,v,s)je( s+<y 2) x / v (p 2sv (-u)dv 


( 1 . 1 ) 


and, for x > a, 


' ,± f o f A l±(- v ) - P ± (-a,v,s) + Fj^X-x^-v^j 
x e~ ( E+CT 1 ) X / V cp ls v ( n)d.v 

± f 1 F 1+ (-x,v,s)e( s - ro l^/ v <» (-^dv. 


( 1 . 2 ) 




vith an equation similar to (l,2) for x < -a. We see then that ve need the coefficients 



... 





0(s) 


v . 2 


x q( z ' s ) "* Of 1 ); 


and, for 0 < (i, v < 1 


<P JSV ( t i) -* 0(s) 


<Pj EV (-d) -*0(1). 


(1*7) 


The Fj+ functions appear with an exponential factor and its behavior 


in the same limit is 


0)/v [f J± (xo,v) +o(|)]to 0 

-*# (l - 8) 


On using Eqs. (1. 7) and (l.8) in Eqs. (1.5) and (l.6) we find that 

^2fc(M"> s ) , . (s*K7i)a/u, /i\ . v 

__--. Pl± (. a ,M,s)e -o(-) (1*9) 


Iq+Cl 1 * s ) 

^ls(^)^ls(M’) 


TFl±(-a,d,6)e- (s+<I l )a /d ± r 2t (a >t i,s)e' {s4<, 2) a /l* 


(I. 10) 


The coefficient A^Cl-O is obtained from the integral equation (1. 3) • 
Since the kernel of this equation is also "the first term of the 

Neumann series solution will give the behavior of Ag+(p) as 
| s | -400. it follows then from Eqs. (1.3) and (1.9) that 
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A2 ±(h) ->F 1± (-a,n,s)e 


-(02-<Ji)s.M _* o m. 


Using Eqs. (1.10) and (i.ll) in Eq. (1.4), we obtain 

[Ait(-n) ± F 1± (-a,|i,s)] ->± F at (a^,s)e“^° 2 “ ai ^ a ^ 

± F 1± (-a,fi,s)e" 2 ( <J 2-< y l) a /^ 


•40 ftV 


These last two results are used in Eqs. (i.l) and (X.2) to get 


n > 0 


±(s +0 2 )e“^ s+0 2) x /^ F 2 f(-x,-n,s) + F 3± (-a,-^,s)e^ 2 “ <Tl ^ a ^, 


±(s + ai)e"( s+cf l) x /^^i±(-x, -|i,s) - Fi+(-a,-n,s) 

± F 2± (a,|i,s)e“^ a 2" cr l^ a ^ 

■ 

± F2+(-a,|i, s)e" 2 ^ a 2" ff l) a /^, (j. > 0 

± (s + ai)e"^ s+a l^ x ^ F 1± (-X,-H,s), M < 0, 


when x > a. For x < -a, has a similar form. Upon using 

♦ 

Eq. (1.8) for the Fj+ functions, we find that Eqs. (1.13) and (l.l4) 
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can be -written as Eqs. (5,l)-(5.4) of the text, where we have used the 
symmetry properties of fj 4 .(x,fx). It can be seen from Eqs. (5.1)-(5*4) 
that is not necessarily 0^. In fact if fj + (xo,|i) 

contains 5(x - xq) then ty s+ (x,u) is 0(1) as |s| -*«, Re(s) = y. 
The parts of ts+( x >^) which are not can 1)6 inverted by 

inspection after a change of variables is made. 

We define |i) for x > a, |x > 0 and all s as 

i us± (*,n) - if* e -(3«l)(x-xo)/^ fl±(x0)Pl)too 

^ ^ -00 

(1-15) 

which gives the portions of Eq. (5*5) which are not 0 for 
| s J ~*a>, Re(s) = y . Now we malice the change of variables 

x - Xq * ut, (I.l6) 

where t > 0 since x > Xq and ji > 0. Equation (1.15) then becomes 
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+ 


+ 


/, 


(x+a)/n 

(x-a)/p. 


noo 

^ (x+a)/ n 


-( s +ai)t „ (it,jji)dt 

e -(sw 2 )t e -(a 2 .a 1 )(a-x)/n f ^ (x . Mt;(Jl)dt 

' (s+ °l )t e-( ff 2-Ol)2a/n ^ . kt)M ) atj 


(I-1Y) 

which is easily seen to he Eq. ( 5 . 6 ) of the text with ^(x^t) given 
hy Eq. (5*9) • For < 0, we use 


x Q - x * |n|t. (l.l 8 ) 

It is seen then that all of the results given as Eqs. ( 5 . 6 ) -( 5 . 10) 
follow. 

Another point to he discussed in this appendix is the contribution 
from the contour Cp (see Eq. ( 5 . 16 ) ) around the right-hand end of the 
branch cut of Vq^(s) as the radius p goes to zero. This branch 
point is located at s = - 0-^(1 - c-jJ so we define 

s + 02,(1 - c x ) = pe 1 ^. (I.I 9 ) 

Here v 01 (s) -»oo as p -»0 as 

V 01 2 £0 ^ ^ -* < * < *• < I - 2 °) 

The branch cut has already been picked so that Vq^(s) is real when 

« 

s is real and greater than -0^(1 - c^). The integral 
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±(x,^)e st ds = p|t s +(x^)e st je^ d<p, (1.21) 


with s given by Eq. (1.19), is zero in the limit p -» 0 if 


lim 

p-»0 



= 0 , 


( 1 . 22 ) 


independent of 9 . As pointed out in section V, the point 
s m -a^(l - ci) may happen to satisfy the eigenvalue condition, 

Eq. (k.6 ) . We assume for the moment that it does not and show later 
what changes are required if it does. The function % s (») -> 0 as 
p -K) as 


°ls(“) p - ^ pei<P ' (1-23) 

so that 

v 01 2q 1s(»)| - (I.*) 

I s = -Oj/l - cj) 

At this "branch point s e A Spj_ so we need to show the behavior of 
all functions given in section 17 as p -»0. This behavior can be 
given in terms of the behavior of v Q1 and 2i s (<») . In the relation- 
ships which follow, quantities which are functions of s will be 
given as 0(v 01 ), 0(l/v Q1 ), 0(fti s ), 0(l), etc, as s -»-cti( 1 - ci). 

For example, 

flls(°°) v 01 2 ** finite ->0(1), (2.22) 
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where we have given the equation number as that from which the 
ship can be seen, 

Xj s (-n) -> 0 ( 1 ), 
xj 3 (^ 02 ) “4 0 (l) , 

♦ 

and 

x js( ±v 0 l) -> 0 (l/v 01 ). 
k s ~> 0 (l). 

h 2 (ai) -> 0(1), ® = ±v 02 , n (0 < \x < 1 ), 

gj(|i) -4 0(1), 
h l(n) “^(^ls) 

and 

h l( ±v 0l) “ > °( 1 / v 0l)* 

S & fo ) -4 0(1) . 

bi± -» 0(v 01 ) . 

B^(-|i) -4 0(1). 
t 2 ±(x,fi,s) -» 0(l) 

and 

tl±(x,n,s) -4 0(v 01 ) . 

Fj±(x,|i,s) -4 0(1), 

P2±(x,±v 02 ,s) 4 0(1) 

and . 

Fl±(x, ^ 01 , s) -4 0( l/v Q1 ) . 





relation 


(A. 8a) 

(c-7) 


(**■•5) 

( 4 . 2 ) 

( 4 . 4 ) 

(4.5) 

( 4 . 1 ) 

( 5 - 10 ) 
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B2±(v) -> 0(l) . 
Bl±(~v) -*0(l). 

a l± ^(l) 

and 

oc^Lj. — » 0(1) . 

Pi± -♦o(l) 

and 

Pat ->o(i). 


(4.10) 

(4.11) 


(4,i4) 


(4.15) 


F±(-a,v 01 ,s) ->0(l/v 01 ) 




F+(-a,v,s) ->0(l). 

(5.10) 

frl± 

tyj am 

- F±(-a,v 01 ,s) J -4 0(v 01 ). 

(4. is) 

&2± 

+ |F 2± (a,v 02 ,s)j ->0(l). 

(4.15) 


fg±( x ^^ s ) ^o(i). 

(4.16) 


tQ+(x,p,s) -*0(v 01 ). 

(4.17) 


We have from these last two relationships and Eq. (1.24) that 

p\|r s +(x,n) -»fp 0(1), (1.25) 


so that Eq. (j.22) is satisfied. Therefore, there is no contribution 
from the integral (1.21) for the case when s * -<Ji(l - cj does not 
satisfy the eigenvalue condition, Eq. (4.6). 
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If the point s » -ax(l - ex) happens to satisfy the eigenvalue 
condition, then the denominator of [ag+ + Fg±(a, v 02 ,s)] , which is 
equivalent to the eigenvalue condition Eq. (4.6), vanishes. It can he 
seen from Eq. (4.15) that the limiting form of this condition ab the 
branch point is ■ 0 and we shall say something about it in the 

last two appendices. If we consider, for such cases, the function 

* s± (x,n) - [&2± + ^ Fg+.(a,v 0 g,s)] t B;fc (x,n), (1.26) 

instead of ^^(x,^) as the integrand of the integral (l.21), then it 
follows that in the limit p -* 0, the contribution from such an integral 
vanishes. The part which has been subtracted from ^ s+ (x,n) in (1.26) 
is considered separately and would appear to have a pole, due to the 
zero in the denominator of [a^j. + F 2± (a,VQ 2 ,s)J . Its contribution 

therefore does not vanish in the limit p -*0; in fact, its contribution 
looks like a discrete residue term. However, the point is not isolated 

(remember that we are considering the branch point of at 

* 

s * -Oh (1 - C-j^) ) so we shall understand that its contribution is 
included in the branch-cut integral term of Eq. (7.1). We shall see 
from the numerical results that such points occur. 
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Bq- Aligns and Procedures for qonnrrtation. of Eige nvalues 
The equations from which tbs time eigenvalue?* 3 & are determined 
€ St* As©? are Eqa. (U,2) and (*.€} of section, It. When 
f% the eorrespcniing equations are Iqs* (*►» t ) 
they -determine what we hate called the pseudo-eigenvalues . Ali or 
these equations are solved numerically using the procedure of references 
1 and 1. As pointed out in section YI, the shove equations can he 
written in terms of the nondimens ional quantities introduced in Iqs. 
(6*1). By making the substitution 


Bfc(n) - 


( v 0i + “) v Pff \ 


n ls^ v 02^ 


is<"> 


) -fe±) 

- (i) V2 6 


(j,i) 


it follows that 3M. (4.2) can toft WWtett for 5 real ana. 

max (—Op + o B , 0) < ? < 1 (thas. -.<s, on that part of the branch out of 

v 0 2 which is not also part of the branch cut of v Q ^) as 


where 


B+(p) - -g(u) 


jft * I v oa| s T + i r 1 ^ j JjL 

— : “ " 2 J n - v ' v + 


U 2 + I v 02 1 S 



** 

All 


UK 


V + 

11 

•X 


g(n) = 


^1^( V 02^ 

%g(*) 


X^(-u) 


u(i - ite 


•2^A/U ^ 


2 . .2 


+ v 


02 




2 2 
c 2 a 2 


(J * i 


«4 


| n i5 (v 02 ) Im J 

fx 2 ^( v 02^ e^ V ° 2 \ 

y o 1£ (») Re i 

[ x l^(v 02 ) V 01 " v 02j 


(J.5) 


and 
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Now we define as 

r>l 

A£± * -2g+ T J B+(v)dv, (j.4) 

which is the eigenvalue condition, Eq. (4.6), if A^ + * 0. 

Equations (J.2) are reduced to two sets (±) of N equations in 
the N unknowns and B.(n i ), 1*1, ,,.,N (see for example, 

ref. 24V, given by 


B ± (di) - -gtn*) 


di6± ± 


'02| 


6, 


N 


L ^i + I v 02 


2 


■ii 


Rj 


J»1 


W. 

rj + 


(j-p) 


where R is the weighting function for the numerical integration 

scheme which is used. Equation (J»4) is written as 

% 

N 

A ?± = -2«± T Y Rj B ± ( Mj ). (J.6) 

Since we must search for values of £ for which A^± = 0, it is seen 
that we must be able to compute all quantities which appear in Eqs. 
(«J.3) for any value of £ in the range (6. 3)* These quantities are 
computed as follows. 
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The functions Vq^ are determined by = 0 an<3, they are 

computed numerically using the Nevton-Raphson Iteration (ref. 2k) on 
the nonlinear equations 


t V 02 | 


tan 


-1 


02 


= s 


(j.9) 


and 


v 0l tanh 


-1 1 


£ + Or 


V 01 <*• 


(J.10) 


R 


The X -functions are computed from the first relationship in Eqs. (A. 8a) 
namely. 


(j.ll) 


X J{( z) - -i- exp 


U. c 1 tn 

nj e (v) 

dv 1 

|2«i J o 

L fi jV v) J 

v - zj 


For £ and z real where z = -p, 0 < n < 1, we have from Eq. (j.ll) 

that 


x gt ( ~ u) 

X 1? (-u) 


= exp 


(i/ 0 l Wt- 




°R'°D' V 1! V 4 : pj 


, (J.12) 


where 


and 


02 (£>v) = tan 


-1 


jtv/2 


©l(k°R' a D' v ) T tan 


-1 


£ - v tanh”- 1 - v 
a R itv/2 


£+Oj) - dpv tanh” 1 v 


(J.13) 


j * • f 
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For £ real and z = v 02 , we calculate the real and imaginary parts 

X 2£^ v 02^ 

of as 


.£( v 02) 


x 2^ v 02^ r l/* v 

Re 7^77 — r = e cos(r 2 /it) 
*l£' V 02' 


X 2t( v 02^ P lA 

Im = e 

X l£^ V 02^ 


sin(r 2 /n), 


where 


r l = " Bl&W)] 2 + ' 


v 02 1 


(j.l4) 


= \ v ° 2 \ f 0 


° E ’ aD ’ V) } v* + |v 0 2 | 2 ’ (J ' 15) 


Integrals in Eqs, (j.12) and ( J.15) are computed as 

1 M 
|e 2 (£/v)- M i, 0%, <*T), v] f ( V )dv B ± |e 2 ( £, V i ) -e^ £,a R ,a D , v-JJffv*), 


(J.16) 


where R i is again the weighting function for the numerical integration 
scheme . 

In all numerical integrations, we used Gauss' Method (ref. 24). 

For integrations in Eqs. (j.5) and ( J.6), the interval (0, l) was 
split into four intervals. 


i 


15 * 

( 0 , 1 ) = (0,0.05) + (0.05,0.1) + ( 0 . 1 , 0 . 9 ) + ( 0 - 9 , 1 . 0 ), (J. 17 ) 

and we used a 10-point Gaussian formula in each subinterval. For 
integrations in Eqs. (j.12) and ( J.15) the interval (O.l) was divided 
as 

(0,1) = (0,0.1) + (0.1,0. 9) + (0.9,0.99) + (0.99 + 0.999) + (0.999,1.0), 

(j.18) 

and in each of these subintervals we also used a 10-point Gaussian 

formula. The subdivision (J.l8) is the same as that used by Kowalska 

(ref. 11) and the X-f unctions calculated here agree with those she 

gives to all figures which she quotes except for the. real and imaginary 

parts of X^(v 02 ). She apparently used r 2 instead of r 2 /jt in Eqs. 

(«J.l4) to obtain the numerical values for the real and imaginary parts 

given in Part II of reference 11. Since her later published critical- 

slab results (ref. 12) agree with those of Mitsis (ref. 20) for a bare 

slab, we expect that this oversight was corrected. 

Conditions (4.7) and (4.8) which determine the pseudo-eigenvalues 

for s e S^ e A S 2 i lead to very similar equations which will not be 

written down. In this region, the real s-axis corresponds to 

0 < £ < -cfj) and such equations need be considered only if -aj) > 0. 

The procedure we use to calculate the eigenvalues is as 

follows. For fixed values of A, cr^ and a^, we select a number of 

^ values in the interval given by (6.3) • For each of these values we 

* 

obtain |Vq 2 | and V 01 from Eqs . (J.9) and (j.10) by iteration 
( Newton -Raphson) . Equations (j.13) are evaluated at each of the 


1 — — 

■. . ..t, 1. • - . W . 
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w x o 1.(^4) 

50 Gaussian integration points, v,. 0 < v. <1. Next, the -=■* 

1 1 

are calculated for each of the 40 Gaussian integration points, 

|ij, 0 < < 1 by using Eqs. (J.l6) in Eq. (j.19). The real and 

v 02 ^ 

imaginary parts of — are computed in the same way from $qs. 

X l^(v02) 


(«J.l4) - (J.l6). Now we can compute g(nj) from Eq. (j.3) at each of 
the 40 points and evaluate all of the coefficients in the two 

1 

sets (±) of N equations in the N unknowns B + ( |Xj ) and 
B_(nj ) (Eqs. ( J.5) ) • These two sets of simultaneous equations are 
solved numerically for B+(nj) which are then used to compute A^ + 
from Eq. (j .6) at the selected values of £. In this way we locate the 
zeros of A^+ approximately. A new set of £ values, located about 
each approximate £ n , is selected and the process is repeated. For 
the present computations, the £ n were located to three figures. 
Discussion of computed results is given in section VI. 

In Appendix G, the thick-slab eigenvalue condition was given as 
Eq. (G.4). We note that g± given by Eq. (j.3) are, within a factor, 
exactly the quantities needed in Eq. (G.4). Therefore the thick-slab 
approximation eigenvalues are obtained from 


8± = 0, / (J.19) 

as would be expected from Eq. ( J.4) . 

The bare slab eigenvalues are obtained when <jp = 0 and it is 
easily shown that in this case Eqs. (J.5) and (J.6) no longer depend 
on <j d * That is, for ctr = 0 these equations do not contain cr D . 
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We have noted in section V and Appendix I that the branch point 
of v Ql located at s = -a^(l - c^) may happen to satisfy the eigen- 
value condition, which can be seen from Eq. (4.15) to be 

a 1± = 0 (J.20) 

when Vq-l -> This point corresponds to £ = -Op + crp and it can be 

shown that Eq. (J.20) then determines values of £ = £ n which do not 
depend on or <jp. That is, if we use £ = -Op + 0 p to eliminate 

CTp from the condition (j.20), then 0p drops out of the equations. 
Equation ( J . 20 ) determines the values of £ at which eigenvalues 
disappear into the right end of the branch cut of Vq^. We also note 
that the limiting form of the pseudo-eigenvalue condition for s = -a^, 
which corresponds to £ = -Op, determines the values of £ where the 
pseudo -eigenvalues disappear into the left end of the branch cut of 
vqi- Such points, as well as those given by Eq. (j.20), are labeled 
with * in Figures 6-9. 
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K. Remarks on Eigenvalue - Branch-Point Coincidence 

In this Appendix we make a few remarks concerning the situation 
when the eigenvalues (or pseudo-eigenvalues) disappear into the branch 
cut of Vq^. This situation is somewhat analogous to that encountered 
by Hintz (ref. 10) for the slab surrounded by pure absorbers. He 
could not say whether a bare -slab eigenvalue (which does not depend 
on o D ) that happened to coincide with -a D belonged to the point 
spectrum or the continuous spectrum for his problem. In the present 
problem, the eigenvalues coincide with a branch point as they disappear 
into the branch cut of Vq^. We have not made a numerical study of the 
branch-cut integral in Eq. (7-1) nor have we looked at the eigenvalue 
condition on another Riemann sheet. We suspect that there is no 
drastic change in the shape of the solution given by Eq. (7*1) when 
an eigenvalue disappears into the branch cut of v Q1 and such studies 
would resolve this point. We pointed out in Appendix J that the 
condition ( J.20), which determines whether or not the point 

t 

s = -aj_(l - c x ), ( £ = -c D + a R ), is a zero of the denominator of 

[ a 2± + | F 2± (a ’ v 02> S 

a R explicitly. In Appendix I, we indicated that the contribution 

from such points should be included in the branch-cut integral since 

it arises from the integration around the branch point. We understand 

then that such a contribution is included in Eq. (7*1) If 

s » -cr-j_(l - c^) happens to satisfy Eq. (j.20). We do not know how such 

zeros of Eq. (j.20) behave or appear in the solution after passing 

♦ 

through the branch point as the material properties are varied* 


) 


given by Eq. (4.13), does not depend on o n or 
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If one considered the problem of a finite slab with symmetric 
reflectors of finite thickness then he might be able to see what is 
happening at the places where the eigenvalues coincide with a oo. 

In such a problem, the solution probably does not inherit the branch 
cut of Vqp but Instead has discrete eigenvalues along it. Even 
though there is another parameter in the problem, the reflector thick- 
ness, one might be able to do a numerical study of all the eigenvalues. 


rnmmmmnmmm 
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